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We present theoretical investigation on the high-frequency collective dynamics in liquids and 
glasses at microscopic length scales and terahertz frequency region based on the mode-coupling 
theory for ideal liquid-glass transition. We focus on recently investigated issues from inelastic-X- 
ray-scattering and computer-simulation studies for dynamic structure factors and longitudinal and 
transversal current spectra: the anomalous dispersion of the high-frequency sound velocity and the 
nature of the low-frequency excitation called the boson peak. It will be discussed how the sound 
mode interferes with other low-lying modes present in the system. Thereby, we provide a systematic 
explanation of the anomalous sound- velocity dispersion in systems - ranging from high temperature 
liquid down to deep inside the glass state - in terms of the contributions from the structural- 
relaxation processes and from vibrational excitations called the anomalous-oscillation peak (AOP). 
A possibility of observing negative dispersion - the decrease of the sound velocity upon increase of 
the wave number - is argued when the sound-velocity dispersion is dominated by the contribution 
from the vibrational dynamics. We also show that the low-frequency excitation, observable in both 
of the glass-state longitudinal and transversal current spectra at the same resonance frequency, is the 
manifestation of the AOP. As a consequence of the presence of the AOP in the transversal current 
spectra, it is predicted that the transversal sound velocity also exhibits the anomalous dispersion. 
These results of the theory are demonstrated for a model of the Lennard- Jones system. 

PACS numbers: 64.70.Pf, 63.50.+X, 61.20.Lc 
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I. INTRODUCTION 



The study of high-frequency collective dynamics in liq- 
uids and glasses at microscopic length scales and ter- 
ahertz frequency region has been a subject of intense 
investigations in the past decade. In particular, recent 
development of the inelastic X-ray scattering (IXS) tech- 
nique has renewed the interest in a long-standing issue of 
sound propagation Q, H Q . It is now generally accepted 
that a well-defined sound-like oscillatory mode - called 
the high-frequency or fast sound - is supported in liquids 
outside the strict hydrodynamic region, down to wave- 
lengths of a few interparticle distances, but with a sound 
velocity larger than the hydrodynamic value. Tradition- 
ally, such increase of the sound velocity upon increase of 
the wave number - called the anomalous or positive dis- 
persion - has been interpreted within the so-called vis- 
coelastic model 0,H|. In this model, relevant memory 
kernel in the Zwanzig-Mori or memory-function equation 
for the dynamics structure factor S q (ui) is modeled by 
an exponential function with a single time constant r 
reflecting structural relaxation in liquids. An increase 
of the sound velocity is predicted at the wave number 
q where the condition w™ ax T = 1 is fulfilled with the 
resonance frequency w™ ax of the dynamic structure fac- 
tor. This condition marks the transition from the low- 
frequency viscous behavior to the high-frequency clastic 
behavior of the liquid, and thus, the positive dispersion 
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effect is ascribed to the transition between a liquid-like 
to a solid-like response of the system. 

According to the viscoelastic model, the absence of 
the positive dispersion is naturally expected in glasses 
where the structural relaxation is basically frozen. How- 
ever, molecular dynamics (MD) simulations have demon- 
strated the presence of the positive dispersion in model 
glasses d 0, H H EH- A first experimental confirma- 
tion has been reported recently from an IXS study on 
vitreous silica [Tj. The observed anomalous dispersion 
implies the presence of some additional "relaxation" pro- 
cess which is active even in glasses, and stimulates an 
extension of the viscoelastic model to fully account for 
the anomalies present both in liquids and glasses. 

Such extension was also motivated by recent IXS stud- 
ies on simple liq uid metals near the melting tempera- 
ture El El El El EH From a detailed line-shape 
analysis of S q (ui), the viscoelastic model was found to 
be unable to account for inelastic peaks in these systems. 
It was found, instead, that a two-time-scale model for the 
memory kernel successfully describes the observed spec- 
tral shapes. Thus, a convincing experimental demonstra- 
tion is provided for the simultaneous presence two re- 
laxation processes with slow and fast characteristic time 
scales, termed structural (labeled a) and microscopic (la- 
beled n) processes. Furthermore, it was pointed out that 
the faster /j, process is the major controller of the pos- 
itive dispersion observed in these systems since it was 
always found that w™ ax T Q > 1 for the slower time scale 
T a , whereas w™ ax r M becomes close to unity at some wave 
number for the faster time scale t m . According to this 
experimental observation, the positive dispersion cannot 
fully be ascribed to the structural relaxation. 
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The nature of these two processes have been studied 
in more detail in Ref. based on the MD simulation 
for a model of lithium. From the analysis of the simu- 
lated S q (uj) with the two-time-scale model, a strong tem- 
perature dependence was observed for the time scale r a , 
corroborating that the slower process is associated with 
the structural relaxation. On the other hand, the faster 
fi process was found to persist both in liquid and glass 
phases with almost no temperature dependence of its re- 
laxation time Ta- The persistence of the p, relaxation 
accounts for the presence of the positive dispersion in 
glasses. Of course, a natural question would be what 
is relaxing at that low T, i.e., what is the microscopic 
nature of the p, dynamics. This was studied in Ref. Q 
based on the computer simulation for a harmonic glass. 
There, it was suggested that the origin of this process can 
be ascribed to the topological disorder, i.e., the decay of 
the memory kernel is due to the dephasing of different 
oscillatory components in the force fluctuations. 

Besides the peak associated with the high-frequency 
sound, the MD [1 [hJ and IXS [H E3 works showed 
the presence of additional low-frequency excitation in 
the dynamic structure factors S q (uj) or in the related 
longitudinal current spectra oc u 2 S q (uj). This second 
low-frequency excitation - also called the boson peak 
- exhibits some general characteristics: it appears in 
the spectra at q values larger than some fraction of the 
structure-factor-peak position, and its resonance posi- 
tion is weakly g-dependent. It was conjectured that 
the low-frequency excitation reflects the transversal dy- 
namics, and the appearance of the transversal mode in 
the longitudinal current spectra was interpreted as being 
due to "mixing" phenomena caused by microscopic dis- 
order |lg|. The main evidence supporting this conjecture 
lies in the simulation results that the low-frequency exci- 
tation appears at the same resonance frequency in both of 
the transversal and longitudinal current spectra, but its 
intensity is more enhanced in the former. Another sup- 
port was inferred from studies on the density dependence. 
The computational work in Ref. |10| demonstrated that 
the resonance frequency of the low-frequency excitation 
in the transversal current spectra increases with increas- 
ing density. Experimentally, it is observed that upon 
densification the boson-peak energy shifts to hi ghe r en- 
ergy and its intensity strongly decreases |2(J, Wi, • 
This parallel of the density dependence was claimed to 
also suggest that the boson-peak-like low-frequency exci- 
tation arises from the transversal dynamics. 

In this paper, we present a microscopic and unified ex- 
planation of the mentioned features of the high-frequency 
dynamics in liquids and glasses - the anomalous disper- 
sion of the high-frequency sound velocity and the na- 
ture of the low-frequency excitation - which have so far 
been investigated somewhat independently. This will be 
done based on the mode-coupling theory (MCT) for ideal 
liquid-glass transition (2^. Originally, the MCT was 
developed to deal with structural-relaxation processes 
which evolve as precursor of the glass transition. On 



the other hand, the glass transition also modifies the 
short-time or the high-frequency dynamics. The study 
of these modifications was the main subject in Ref. |24| . 
It was shown there that the strong interaction between 
density fluctuations at microscopic length scales and the 
arrested glass structure causes an anomalous-oscillation 
peak (AOP), which exhibits the properties of the bo- 
son peak. As will be demonstrated in the present study, 
the AOP persists in liquid states as well as in glass 
states. It will be discussed how the sound mode in- 
terferes with other low-lying modes present in the sys- 
tem. Thereby, we provide a systematic explanation of 
the anomalous sound-velocity dispersion in systems - 
ranging from high temperature liquid down to deep in- 
side the glass state - in terms of the contributions from 
the structural-relaxation processes and from the vibra- 
tional modes building the AOP. A possibility of observing 
negative dispersion - the decrease of the sound velocity 
upon increase of the wave number - is argued when the 
sound-velocity dispersion is dominated by the contribu- 
tion from the vibrational dynamics. We also show that 
the low-frequency excitation, observable in both of the 
glass-state longitudinal and transversal current spectra 
at the same resonance frequency, is the manifestation of 
the AOP. This seems to be an alternative interpretation 
to the mixing phenomena mentioned above. As a natural 
consequence of the presence of the AOP in the transver- 
sal current spectra, it is predicted that the transversal 
sound velocity also exhibits the anomalous dispersion. 

The paper is organized as follows. In Sec. [HJ the ba- 
sic MCT equations of motion are formulated, and the 
model for the demonstration of our theoretical results is 
specified. Section II I II deals with the anomalous sound- 
velocity dispersion of the high-frequency sound, while in 
Sec. |W] we discuss features of glass-state longitudinal and 
transversal current spectra at low frequencies. Section Ivl 
concludes the paper. 

II. BASIC THEORY 

A. MCT equations of motion 

The basic quantity describing the equilibrium struc- 
ture of a simple system is the static structure factor 
S q = (| Pq-\ 2 } defined in terms of the density fluctuations 

for wave vector q, pq = J2iLi exp(i<f • r*j )/VN. Here, N 
is the total number of particles in the system distributed 
with the average density p, and f\ denotes the position 
of the ith particle. For homogeneous and isotropic sys- 
tem, which we assume throughout the paper, the static 
structure factor depends only on the modulus q = \q\. 
The most relevant variables characterizing the structural 
changes as a function of time t are the density correlators 
4> q {t) = (pq-(t)* Pq(0)} / S q which are normalized to unity 
at t = 0. The short-time asymptote of these functions is 
given by cj> q (t) — 1— (1/2)Q^ 2 +- ••, where fl 2 denotes the 
square of the characteristic frequency Q 2 = q 2 v 2 /S q 
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Here v = ^k^T/M with Boltzmann's constant fee de- 
notes the thermal velocity of particle of mass M at tem- 
perature T. In the small wave-vector limit, one obtains 
= v o c l + 0(q 3 ) with the isothermal sound velocity 
vq = v I -J S q =o ■ Within the Zwanzig-Mori formalism [f| 
one can derive the following exact equation of motion 

afo,(f)+n^ g (f)+fi* f dt'm q {t-t')d tl <t> q {t') = 0, 
Jo 

(la) 

in which the relaxation kernel m q (t) is a correla- 
tion function of fluctuating forces. Let us introduce 
Fourier-Laplace transformations to map <j> q , m q , and 
similar functions from the time domain onto the fre- 
quency domain according to the convention 4> q (uj) = 
i J Q dt exp(iu)t)(j> q (t) = (f>' q (oj) + i(f>g(u). Equation fla|) 
is equivalent to the representation 

= -l / { w - n* / [u + n 2 q m q (uj) } }. (lb) 

Within the MCT, the kernel m q (t) is given by a mode- 
coupling functional of density correlators describing the 
cage effect of dense systems, 

m q (t) = F q [<t>{t)}. (2a) 

Here the mode-coupling functional reads F q [f] = 
1/(2tt) 3 / dkV(q;k,p) f k fp with p = q-k. The cou- 
pling coefficients V(q;k,p) are determined by the equi- 
librium structure. Within the convolution approximation 
for triple correlations, they are given by V(q;k,p) = 
pS q SkSp [q ■ (kck + pc p )] 2 I (2g 4 ), where c q denotes the 
direct correlation function related to S q via the Ornstein- 
Zernike equation, pc q = 1 — 1/S q 0. For details, the 
reader is referred to Ref. |2^|. In the following, wave- 
vector integrals will be approximated by Riemann sums. 
The mode-coupling functional can then be written as a 
quadratic polynomial 

k.p 

where the indices run over the discretized wave-number 
grids. The coefficients V q ^ P are trivially related to 
V(q ;k,p). The q — > limit of the functional reads 

k 

with V k = {pS /An 2 )k A S 2 [ C 2 + 2(kc' k )c k /3 + (kc' k ) 2 /5}. 

Equations Q and (0) are closed, provided equilibrium 
static quantities are known as input. In the following, the 
results will be demonstrated for the Lennard- Jones (LJ) 
system: the interaction between two particles is given by 
the LJ interaction V(r) = 4e L j{(cr L j/r) 12 - (ct L jA) 6 }. 
The first MCT work for the LJ system has been done in 
Refs. pM l26j | with the optimized random-phase approx- 
imation for S q In the present study, S q is evaluated 



within the Percus-Yevick approximation y. Wave num- 
bers will be considered up to a cutoff value q* — 80/ctlj, 
and they are discretized to 400 grid points. From here on, 
all quantities are expressed in reduced units with the unit 
of length ctlj, the unit of energy 6lj (setting fee = 1), and 
the unit of time (Mctlj/clj) 1 ^ 2 - The dynamics by vary- 
ing T shall be considered for a fixed density p = 1.093, 
except for Sec. IIV 61 where the dynamics at a lower den- 
sity p = 0.95 is discussed. (We notice that the triple point 
of the LJ system is at p* » 0.85 and T* w 0.68 |27|.) 

B. Ideal glass states 

The specified model exhibits a fold bifurcation [23| . 
For small coupling constants V q ^p, the correlators <j) q (t) 
and the kernels m q {t) decay to zero for long times, <fr q (t — > 
00) = and m q (t — > 00) = 0. In this case, the spectra 
4>q(uj) and rtiqiw) are continuous in ui. In particular, 
there hold lim^^o ^4>q( UJ ) — and lim^^o u>m q (ui) = 0. 
Density fluctuations created at time t = disappear 
for long times and the same holds for the force fluctu- 
ations as expected for an ergodic liquid. For large cou- 
pling constants, on the other hand, there is arrest of 
density fluctuations for long times: <p q {t — > 00) = f q , 
< f q < 1. Thus, nonergodic dynamics is obtained 
in which the perturbed system does not return to the 
equilibrium state. Similarly, there is arrest of the force 
fluctuations, m q (t — * 00) = C q > 0. The nonergodicity 
parameters f q and C q are connected via |28| 

f q = C q /(l + C q ), C q =T q [.f]. (3) 

For this strong-coupling solution, the kernel exhibits a 
zero- frequency pole, lim^^owwigfw) = — C q . The fluc- 
tuation spectrum 4> q (u>) oc S q (u>) exhibits a strictly elas- 
tic peak: 4> q (u>) — nf q 8(oj) + regular terms. This is the 
signature for a solid with f q denoting its Debye- Waller 
factor. Hence, the strong coupling solution deals with a 
disordered solid; it is a model for an ideal glass state. If 
one increases the coupling constants smoothly from small 
to large values, one finds a singular change of the solu- 
tion from the ergodic liquid to the nonergodic glass state, 
i.e., an idealized liquid-glass transition. For simple- liquid 
models, the transition occurs upon cooling at some criti- 
cal temperature T c or upon compression at some critical 
density p c [2^. For the LJ model under study where 
the temperature is considered as a control parameter, 
one finds T c « 1.637 for p — 1.093 within the Percus- 
Yevick approximation for S q . If T decreases below T c , 
f q increases above its value at the critical point, called 
the critical nonergodicity parameter or the plateau f q , as 
demonstrated in Fig. ^ 

C. MCT liquid-glass-transition dynamics 

As precursor of the ideal liquid-glass transition, MCT 
predicts the evolution of the glassy dynamics which are 
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FIG. 1: Debye- Waller factor f q and one-tenth of the static 
structure factor S q at p = 1.093 for T = 0.5 (solid lines), 
T = 0.8 (dashed lines), and T = T c » 1.637 (dotted lines). 



stretched over many decades. Such MCT scenario has 
been comprehensively discussed for the hard-sphere sys- 
tem in Refs. [2|||3(]]. This subsection compiles the MCT 
universal predictions for the dynamics near the liquid- 
glass transition to an extent which is necessary for un- 
derstanding the present article. 

We start from introducing some concepts to describe 
the MCT-liquid-glass-transition dynamics |2^|. In the 
space of control parameters, a smooth function a is de- 
fined near the transition point, called the separation pa- 
rameter. Glass states are characterized by a > 0, liquid 
states by a < 0, and a = defines the transition point. 
When only the temperature T is considered as a control 
parameter near the transition point, one can write for 
small distance parameters e = (T c — T)/T c : a = Cye, 
Ct > 0. In addition, the transition point is character- 
ized by a time scale to an d by a number A, < A < 1. 
The scale to specifies properties of the short-time tran- 
sient dynamics, and A is called the exponent parameter. 
The latter determines a certain number B > 0, the crit- 
ical exponent a, < a < 1/2, and the von-Schweidler 
exponent b, < b < 1. For the specified LJ model at 
p = 1.093, we have C T = 0.226, A = 0.697, a = 0.328, 
b = 0.646, B = 0.669, and t = 0.0125. 

Let us consider the correlator 4>x{t) = (X(t)*X(0)) 
of some variable X coupling to the density fluctuations. 
Its nonergodicity parameter fx — <px(t ~* °°) exhibits a 
square-root singularity near the transition, 



fx =f x + W<V(1-A), a > 0, 



0. 



(4) 



The critical nonergodicity parameter f x and the critical 
amplitude hx arc equilibrium quantities to be calculated 
from the relevant mode-coupling functionals at the criti- 
cal point T c . At the transition, the correlator exhibits a 
power-law decay that is specified by the critical exponent 
a. In a leading-order expansion in (l/i) a , one gets 



For small values of e, there is a large time interval, 
where <fix(t) is close to its critical nonergodicity param- 
eter f x . Solving the equations of motion asymptotically 
for this plateau regime, one gets in leading order in the 
small quantities 4>x{t) — fx the factorization theorem: 



4>x(t) = f x + h x G(t). 



(6a) 



The function G(t), called the (3 correlator, is the same for 
all variables X, and describes the complete dependence 
on time and control parameter via the first scaling law 



G(t) 



'g±(t/t a ), a^0. 



(6b) 



Here, t a = t a /\a\ i - ' ' is the first critical time scale. 
The master functions g± (i) are determined by A |3l| . For 
t <C 1, the master functions are the same on both sides of 
the transition point, and one gets the power law g±(t — > 
0) = l/t a so that Eq. J^J is reproduced for fixed large t 
and a tending to zero. For f > 1, one gets for glasses 
g + (i > 1) = 1/Vl — A, so that Eq. Q is reproduced. On 
the other hand, for liquids, one finds for the large rescaled 
time t: g_(i -»oo) = -Bt b + 0{l/i b ). Substituting this 
result into Eqs. JBJl, one obtains von Schweidler's law for 
the decay of the liquid correlator below the plateau f x : 



4>x{t) = fx-hxit/t'J, t ff «t, 



-0. (7) 



The control-parameter dependence is given via the sec- 
ond critical time scale t' a = to B -1 / 6 / \a\^ with 7 = 
(l/2a) + (1/26). The dynamical process described by 
the cited results is called the /3-process. The (3 correlator 
G(t) describes in leading order the decay of the correlator 
towards the plateau f x within the interval to ^it^it a . 
For t ^> t a , the glass correlator arrests at fx, while the 
liquid correlator exhibits von Schweidler's law. 

The decay of the liquid correlator below the plateau f x 
is called the a-process. For this process, there holds in 
leading order for a — > — the second scaling law, 4>x(t) = 
<px (t) with t — t/t'p, which is also referred to as the super- 
position principle. The control-parameter-independent 
shape function <fix(t) is to be evaluated from the mode- 
coupling functionals at the critical point. For short 
rescaled times t, one gets <j>x(t) = fx — hxi + 0(t 2b ), so 
that Eq. (J2J is reproduced. The ranges of applicability of 
the first and the second scaling laws overlap; both scaling 
laws yield von Schweidler's law for t a <C t <C t' a . 



D. Description of the glass-state dynamics 

For later convenience, we introduce a reformulation of 
the MCT equations of motion which is more appropriate 
in handling the dynamics in glass states 23]. We map 
the density correlators 4> q to new ones <p q by 



ox[ 



(t)=f x +h x (to/tr, 



0, (t/to) 



00. 



(5) 



<j> g (t) = f q + (i-f q )<f>M 



(8a) 
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This amounts to dealing only with the decay relative to 
the frozen amorphous structure described by f q . Intro- 
ducing new characteristic frequencies Q, q by 

n» = n»/(i-/ g ), (8b) 

one obtains for short times <f> q (t) = 1 — (l/2)tl 2 t 2 + ■ ■ ■ 
in analogy to the one for <fr q (t). In the small wave- vector 
limit, there holds Cl q = ioq+0(q 3 ) with vo = Vo/y/l — fo- 
The modification of the sound velocity reflects the re- 
duced compressibility in nonergodic glass states. Substi- 
tution of these results into Eq. (|lafl reproduces the MCT 
equations of motion with cf> q , Q q , and m q replaced by cf> q , 
Q q , and rh q , respectively. Here, the new relaxation kernel 
is related to the original one by 

m q (t)=C q + (l + C q )m q (t). (8c) 

The new correlator has a vanishing long-time limit, 
\im t oc 4> q (t) = 0, and correspondingly the Fourier- 
Laplace transform exhibits a regular zero-frequency be- 
havior lim w _>o uj<f> q (uj) — 0. Combining Eqs. and 
||ScJ, one also concludes that lim^oo fh q (t) = and 
lim^^o w?fig (u) = 0. Since the equation of motion does 
not change its form, one gets as the analog of Eq. (|lb|) 

$ q (oj) = -1/{co-C1 2 /[lo + & q m q (uj) } }. (8d) 

The reformulated equation of motion for <p q (t) can be 
closed by finding an expression for the new kernel rh q (t) 

as a new mode-coupling functional T q . One finds com- 
bining Eqs. ijgal) and JSc} with Eqs. l|2al) and lj2b"|l that 
T q is given by a sum of linear and quadratic terms 

m q (t) = =F (1) W)] +& 2) [*(t)], (9a) 

where 

^[/]=EV^ ^'[/VEW^. (9b) 

k k,p 

with renormalized coefficients 

V q , k = 2(1- f q )J2v q ,k P (l-fk)f P , (9c) 
p 

V g ,k P = (l-f q )V q>kp (l-f k )(l-f p ). (9d) 

Correspondingly, one finds from Eq. (Pc|) for the zero 
wave-number limit of the kernel m(t) — m q =o(t) 

m(t) = m (1) (i) +m (2) (f), (10a) 
in which linear and quadratic terms are given by 

k k 

(10b) 



with new coefficients 

Vt 1] = 2(l-f )V k f k (l-fk), (10c) 
Vk 2) = (1- /o)V fc (l- fk f. (10d) 

As a result, equations of motion are produced, which 
are of the same form as the original MCT equations 
and J2J) • But in addition to the quadratic mode-coupling 
term, there appears a linear term, and the original mode- 
coupling coefficients are renormalized to hatted ones. 
The linear term describes interactions of density fluctu- 
ations with the arrested amorphous structure, while the 
quadratic one deals with two-mode decay processes. 

E. Longitudinal and transversal current correlators 

Current correlators are also relevant variables de- 
scribing the dynamics of dense systems, in particular, 
when one is interested in vibrational properties. These 
are defined in terms of the current fluctuations, j'2 = 

J2iLi v i ex p(*9 • ?i )/y/~N, where a (= x, y, or z) refers to 
the spatial component and vf denotes the velocity of the 
ith particle. Since (v™Vj) = 5ijS a /3 v 2 , one gets the static 
correlation = 5 a p v 2 . The current correlators for 

isotropic system can be represented by two independent 
functions which depend only on the modulus q: 

oi«*jf (o)>/« 2 = rf<i>\{t)+[5 a p-ef]<t> T q {t). (ii) 

Here q a denotes the a component of the unit vector along 
q. The functions 4> q (t) and <j)J(t) are called the longitudi- 
nal and transversal current correlators, respectively, and 
they are normalized to unity at t = 0. 

1. Longitudinal current correlator 

The density fluctuations and the longitudinal cur- 
rent fluctuations are related via the continuity equation, 
d t pq- = iq ■ js; and there holds d 2 <j) q (t) = —fl 2 (f> q (t). 
Their Fourier-Laplace transforms therefore satisfy ui\l + 
u>(j) q (uj)] — Q 2 (f>^(ui), and one obtains from Eq. (|lb|) for 

the longitudinal current spectrum 0^ (cu) 
uj 2 n 2 m"(uj) 

^ 1 '~ [w 2 -02A g ( W )] 2 + [o;02 m ^(a;)] 2 ' 1 ' 

with A 9 (cj) = 1 — u) m' q (uj). 

The continuity equation holds also in glass states, 
d 2 4> q (t) = — f] 2 , (j) ^(t) in terms hatted variables defined 
in Eqs. (|8a|) and l|8bl> . and their spectra are related via 

u> 2 4>' q , (cj)=£l 2 q $"(u). (13) 

We therefore have the same expression for </>^ (ui) for 
glass states as in Eq. <|12[) but with all quantities in the 
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right-hand side replaced by hatted ones: 



(glass), 
(14) 

with A q = l — ujm' q (uj). Thus, the same relaxation kernel 
m q or rh q as the one for the density correlator describes 
the dynamics of the longitudinal current. 



2. Transversal current correlator 

The MCT equations of motion for the transversal cur- 
rent correlator 4> q (t) have been derived in Ref. [2ij. The 
Zwanzig-Mori equation is given by 

^(t) + (n^) 2 f dt'm q ? {t-t')<f r q {t')=^ (15a) 



(15b) 



and its Fourier-Laplace transform reads 



1 



c + (^)2mJ( W )' 



with (fiT) 2 = q 2 l 



The relaxation kernel mT (t) is a cor- 



relation function of transversal fluctuating forces. Within 
the MCT, the transversal fluctuating forces are approxi- 
mated by their projection onto the subspace spanned by 
the density products, and with the factorization approxi- 
mation the kernel is given by a mode-coupling functional 
of the density correlators 



k,p 



Vq,kp fkfv 



(16a) 



1 1 (2q 2 ) with 



Here Vj kp = pS k S p [e T ■ (kc k + pc p ) ] 
denoting a unit vector orthogonal to q. The q — > limit 
of the functional is given by 



(16b) 



zero for long times, m q (t — > oo) = C q > 0. We therefore 
write in analogy to Eq. l(5c}l 



(17) 



and introduce m^(t) for the description of the glass-state 
transversal current dynamics. The new kernel satisfies 
lim^oo (t) = and lim^^o ujfn^ (uj) = 0. Substi- 
tuting the Fourier-Laplace transform of Eq. (|17fl into 
Eq. I|15b(l yields for the transversal current spectrum 



Cm 



(18) 



where f2j = J , 



T' 



' and A q = 1 — u> m q 



The reformulated equation of motion for the glass-state 
transversal current dynamics can be closed by finding 
an expression for m^(t) as a new mode-coupling func- 

One finds from Eqs. l|5a)l. I|16a[) . and l(T7|) 



tional JFj 
that.?? = A T(1) 



respondingly, there holds ?71t(£) 



T q in analogy to Eq. Ij9a(l . Cor- 

(£) for 



;/C '(t) 



rn 



(2) 



TOt(*) = ™J = o (0- Explicit expressions for the new mode- 
coupling functionals shall be omitted here for brevity: 
they take the same form as in Eqs. JJJJ and 1)10(1 . but 
with coupling coefficients replaced by the ones for the 
transversal current correlators. 

The mathematical form of Eq. I|18l) for glass states 
is identical to that for the longitudinal current given in 
Eq. ifTljl. In particular, </) q (uj) for glass states vanishes 
proportional to uj 2 for small frequencies in contrast to the 
one in liquid states. Thus, the transversal current cor- 
relators exhibit a drastic difference between liquid and 
glass states [22|: the long wave length correlators in liq- 
uids describe diffusive processes, while the ones in the 
glass states describe the propagation of transversal sound 
waves as expected for an isotropic elastic continuum. 



3. Velocity correlator 

Finally, we briefly summarize here for later discussion 
the MCT equations of motion for the normalized velocity 
correlator $>(t) = (v s (t) ■ v s (0)) / (3v 2 ) defined with the 
velocity v s of a tagged particle (labeled s). An exact 
equation for this quantity is given by 



d t ^(t) + v 2 dt'm s (t-t')y(t')=0, (19) 
Jo 

and the MCT expression for the kernel m s (t) reads 



with V? = (p/60tt 2 )/c 4 [c' k S k ] 2 . These MCT equations of 
motion for qb^ (t) can be solved with the knowledge of the 
static quantities and of the density correlators. 

In glass states, mj(uj) acquires a zero- frequency pole 
since the transversal force fluctuations do not decay to with V s 



m s (t) 



v k s ^ s k (t)Mt), 



(20) 



(p/6n 2 )k*c 2 S k [2 



Here 



;(*) 



(p|(t)*p|.(0)) with pt = exp(iq ■ f s ) denotes tagged- 
particle density correlator. The Zwanzig-Mori equation 
for <j> q (t) has the same form as Eq. I)la[l with 



and til replaced by 



and (fl 



s\2 



qi m q-> 

q 2 v 2 , respec- 



tively, and the MCT kernel is given by the functional 
m s q (t) = Ek^lkM^Mt) with V q * kp determined by 
S q [HEl- For glass states, the MCT equations for tf(t) 
can be reformulated in terms of the hatted kernel rh s (t) 
defined in analogy to Eq. (|17fl via 



m s (t) = C s + C s m s (t), 



C s = lim m s (t). 

t — >oo 



(21) 



7 



60 
40 
20 


60 
40 
20 



~i 1 1 r~ 



{a)T=10 / ooc^^c 




12 3 4 5 
4 1 1 — ^-l 



(b) T=4.0 



oocooo^ >6 



o°° "SO 



(c)T=1.8 



10 



1 1 1 f— 

- - ■ V q 

oO°°°o 



12 3 4 5 
1 , 




FIG. 2: Peak positions 



(circles) of the longitudi- 



nal current spectra cj>^ (u>) as a function of q at p — 1.093 
for temperatures indicated in each panel referring to liquid 
states. Dashed lines denote the hydrodynamic dispersion law 
qvo with the isothermal sound velocity vq- Dash-dotted lines 
show the linear dispersion law qvoo with the infinite-frequency 
sound velocity v x . The insets exhibit apparent sound veloci- 
ties v q — Dq ma:x /q (circles) as a function of q along with hor- 
izontal lines denoting Vo (dashed lines) and Voo (dash-dotted 
lines). Dotted lines in (c) refer to the linear dispersion law 
qvZ> and the velocity discussed in the text (cf. Sec. llIIDl , 
Long-dashed lines through circles in (c) are from the solution 
to Eq. (I23H . while solid lines are based on Eqs. 11248 . 



III. ANOMALOUS DISPERSION OF 
SOUND- VELOCITY 

A. Dispersion relation 

Circles in Fig. [21 show the dispersion relation - the 
peak positions iv^ max of the longitudinal current spectra 

cj)^ (lu) as a function of wave number q - for three repre- 
sentative temperatures referring to liquid states. (We no- 
tice that, in experimental and computer-simulation stud- 
ies for the dispersion relation cited in Sec. U it is the 



peak positions of the longitudinal current spectra that 
are usually reported. The peak positions of the dynamic 
structure factors oc 4>'g(u>) and of the longitudinal cur- 



rent spectra (fit (iS) nearly coincide as can be inferred, 
e.g., from Fig. 1121 In the present article, we are primar- 
ily interested in the wave-number regime q < 4 where 
w Lmax i ncreases with q. In this pseudo first Brillouin 

zone, the peaks ui^ max are associated with the collective 
mode.) The dashed line in each panel denotes the hydro- 
dynamic dispersion law qvo with the isothermal sound 
velocity vq |33|. whereas the dash-dotted line exhibits 
the linear dispersion law qv^ defined in terms of the 
infinite- frequency sound velocity v^ (see below). The 
inset in each panel shows an apparent sound velocity 
v q = ui q jmax /q (circles) as a function of q, along with 
horizontal lines denoting vq (dashed line) and (dash- 
dotted line). The increase of the sound velocity v q from 
vo towards Voo upon increase of the wave number g, which 
is observable for q < 2 in all the insets of Fig. [21 is called 
the anomalous or positive dispersion. 

There are subtle variations in the anomalies depend- 
ing on the temperature. At high temperature T = 10, 
the whole increase of the sound velocity is observable in 
Fig. [2{a) starting from vq at small q up to close to v^ 
at q 2. This feature is altered if the temperature is 
decreased. For T = 4.0 (Fig. Gib)), the sound velocity 
v q for small wave numbers is still located above the hy- 
drodynamic value vo. For T = 1.8 (Fig.[2Ic)), this "gap" 
becomes larger, and the sound velocity v q in the small- 
q limit seems to approach another limit denoted by the 
dotted line which will be defined below. 

Corresponding results for representative glass states 
are shown in Fig. [31 But here, the hydrodynamic dis- 
persion law qvo is determined by vo introduced in con- 
nection with Eq. I|8bjl (see also below). It is clear from 
Fig. [3 that the anomalous sound-velocity dispersion is 
present also in glass states. Again, there are subtle vari- 
ations in the appearance of the anomalies. For T = 1.47 
(Fig-E^b)), the sound velocity v q approaches the hydro- 
dynamic value vo rather rapidly for q — > 0. On the other 
hand, for T = T c (Fig. 01a)), the approach towards £>o is 
retarded. For the deep-in-glass state T = 0.5 (Fig-EI c )); 
a new feature shows up, and "negative" dispersion is ob- 
servable for q < 1, where the sound velocity v q is smaller 
than the hydrodynamic value wo- 

The results shown in Figs.[2Ib),[2Ic), and^a) indicate 
that the increase of the sound velocity from the hydrody- 
namic value has already started at smaller wave-number 
regime which cannot be adequately resolved with the lin- 
ear q axis. In addition, the presence of the negative dis- 
persion shown in Fig. |2Ic) implies that the anomalous 
sound-velocity dispersion cannot be accounted for solely 
in terms of relaxation processes since those processes 
would lead only to the positive dispersion as known, e.g., 
from the viscoelastic model. 

In the following, a microscopic and unified understand- 
ing of the mentioned anomalies shall be attempted. This 
will be done based on the generalized hydrodynamic de- 
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FIG. 3: The same as in Fig. [5] but here temperatures re- 
fer to glass states, and the hydrodynamic dispersion law qvo 
is determined by vo defined in the text. Long-dashed lines 
through circles in (a) are from the solution to Eq. 12311 . while 
solid lines are based on Eqs. EH - 



scription [24J. In this description, the memory kernel 
m q (u>) is approximated by its q — > limit, while the full 
u> dependence of m{uj) — m q= o(u/) is retained. In this 
manner, one gets from Eq. Q12[l [34J 



to 2 Q 2 q m"{u) 



[lo 2 - n 2 A(cj)] 2 + [u n 2 m"{u>)] 2 ' 



(22) 



with A (a;) = 1 — u)m'(u>). In lowest order, the resonance 
frequency uj q max can be obtained as a solution to 



,L max 



(23) 



The long-dashed lines through circles in Figs. and 
EUa) show the dispersion law and the sound velocity based 
on the solution to this equation. One understands that 
the lowest-order solution describes the results from the 
full MCT solutions (circles) fairly well for the whole q- 
range shown in the figure. Since we are interested in the 



small wave-number region q < 2 where the anomalous 
sound- velocity dispersion is present, a further simplifica- 
tion is possible: one can replace Vt q by its leading-order 
contribution, Q q ss qvQ. This yields 



,L max 



-qv Q ^/A(^w, L ™), (24a) 
which can be rewritten for the the sound velocity v q as 



Vo 



^A(uj = qv q ). 



(24b) 



The solid lines in Figs.^c) and 01a) exhibit the disper- 
sion law and the sound velocity based on the solutions to 
these equations. One understands that no serious error is 
introduced with the use of Eqs. (|24l) as far as the regime 
q < 2 of interest is concerned. 

Before proceeding, let us show that Eqs. (|24|l reproduce 
known results [U, U in the limit of small and large reso- 
nance frequencies. The hydrodynamic sound is obtained 
by taking the oj — > limit of A(uj) in Eqs. (J23J. One ob- 
tains the dispersion law uj q max = qvoy^l + m(t — > oo) 
since lim^^o U}m ' '(<*>) == — li m t— ►<» m(t). For liquids 
where the memory kernel relaxes to zero for long times, 
m(t — > oo) = 0, the dispersion law w^ max = qvp is ob- 
tained with the isothermal sound velocity v 33]. For 
glasses, on the other hand, the kernel does not decay to 
zero for long times, m(t — > oo) = C, leading to the dis- 
persion law w^ max = qvo with vo = UoVl + C. (This 
agrees with vq = Vo/y/l — Jo given below Eq. I|8b(l be- 
cause of Eq. (jjJJ.) These results have been included with 
dashed lines in Figs. [21 and [3] In the high-frequency limit, 
one obtains, using lim^^oo A(w) = 1 + lim t _>o m(i), the 
dispersion law max = qv^ with the infinite- frequency 

sound velocity Voo = v y/l + m(t = 0). This holds for 
both liquid and glass states, and has been included with 
dash-dotted lines in Figs. |2] and 

Thus, on the basis of Eq. I|24b|) . it is clear that 
the anomalous sound-velocity dispersion from the low- 
frequency hydrodynamic value (vo for liquids and Vo for 
glasses) towards the infinite-frequency one upon in- 
crease of the wave number q is controlled by the de- 
tailed frequency dependence of A(cj). This in turn is 
determined by the relaxation of the memory kernel m(t). 
Within the MCT, m(t) is given in terms of the mode- 
coupling functional as m(t) = ^[^(t)] [cf. Eq. ij2*c|) ]. In 
the following two subsections, the details of the relax- 
ation of m(t) and of the frequency dependence of A(w) 
will be investigated. These results will be employed in 
Sec. IIII Dl for a systematic study of the anomalies. 



B. Structural-relaxation processes 

Figure 0] exhibits the evolution of the dynamics of 
m(t) upon decrease of the temperature T. The curve for 
T = 10, which is quite higher than T c w 1.637, decays 
rapidly and nearly exponentially to zero. By lowering the 
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FIG. 4: The memory kernel at q — 0, m(t) = m q =o(t), as a 
function of log 10 t at p = 1.093 for temperatures indicated in 
the figure. Curves with labels x = 2 and 3 are calculated for 
distance parameters e = (T c — T)/T c = =pl0 _a: for liquids (e < 
0) and for glasses (e > 0). Dotted line with label c denotes 
the kernel at T c . Horizontal line marks the plateau height C c . 
Dashed lines show the MCT asymptotic-law results based on 
Eq. @3 for T = 1.8 and x = 3 (e < 0). The kernels for 
T = 0.5, 0.6, and 0.8 refer to the right vertical axis. 



temperature, the dynamics becomes slower and stretched 
due to the development of the cage effect. For temper- 
atures close to but above T c , the kernel m(t) exhibits a 
two-step relaxation: the dynamics towards the plateau 
C c , followed by the final relaxation from the plateau to 
zero. The dynamics that occurs near the plateau is re- 
ferred to as the process, whereas the final relaxation 
is called the a process (cf. Sec. I11C|) . The (3 process 
for small \ T — T c \ can be well described by the MCT 
asymptotic formula (|(jafl specialized to m(t) 



m(t) = C c + DG{t), 



(25) 



as exemplified with the dashed lines in Fig. (C c = 
0.982 and D = 1.925 for the model under study.) At the 
critical point T c , the kernel approaches the plateau in 
a stretched manner, as shown by the dotted curve with 
label c in Fig. This process, called the critical decay, 
is described by a power law [cf. Eq. (J5J] 



m(t) =C C + D(t Q /t) a . 



(26) 



Decreasing T below the critical temperature T c , the val- 
ues for the long-time limits C increase. These limits are 
approached exponentially fast for T^T C 35]. If T is far 
below T c , one observes enhanced oscillatory features for 
times near log 10 t ~ —1.2 as can be seen from the curves 




-2-10 1 



FIG. 5: Susceptibility spectra um"(ai) for the memory ker- 
nels m(t) shown in Fig. [I] for liquid states (a) and for glass 
states (b). The insets exhibit the double logarithmic repre- 
sentation of the spectra for T close to T c . The meaning of 
the labels x = 2 and x — 3 is the same as in Fig. 2] Dotted 
line with label c denotes the susceptibility spectrum at T c . 
Dashed line in the upper inset shows the MCT asymptotic- 
law result based on Eq. 1271 for x = 3 (e < 0). Dash-dotted 
lines in the insets exhibit the hydrodynamic linear law ~ u). 



for T = 0.8, 0.6, and 0.5 in Fig. H Such vibrational 
dynamics reflects the anomalous-oscillation peak (AOP) 
discussed in Ref. [24j]. It is caused by strong interaction 
between density fluctuations at microscopic length scales 
and the arrested glass structure. Features of the AOP 
will be discussed in the next subsection. 

Figure [51 shows the susceptibility spectra ujm"{uf) of 
the kernels shown in Fig. The spectrum at T = 10 
exhibits only one peak located at high frequencies, to be 
called a microscopic peak. By lowering T, an additional 
peak emerges at low frequencies reflecting the evolution 
of the structural-relaxation processes. The lower fre- 
quency peak, whose position decreases strongly by low- 
ering T, is due to the a process, and is called the a peak. 
It is separated from the microscopic peak by a minimum, 
called the j3 minimum. The appearance of the minimum 
is caused by the crossover from a power law ~ t~ a to 
another one ~ —t b which occurs in the f3 regime (cf. 
Sec. IllCf) . For small \ T — T c \, the /3 minimum is well 
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described by the formula which follows from Eq. (|25|l , 

um"(w)=DwG"(u), (27) 

as exhibited by the dashed line in the inset of Fig.|SJa). 
At T — T c , the critical decay l|26(l leads to a sublinear 
susceptibility variation 

cum"(uj) = D sin(7ra/2)r(l - a)(ujt ) a , (28) 

as shown by the dotted line with label c. Here, T(x) 
denotes the Gamma function. Below T c , the final expo- 
nential decay towards C |35| implies the hydrodynamic 
law in the small-w limit, Lum"(cu — > 0) oc lu. But, near 
T c , there is a frequency interval where the spectrum is 
described by the sublinear critical spectrum (|28|l . This 
produces a "knee" at some frequency where the crossover 
from the linear low-frequency spectrum u;m"(uj) oc a; to 
the sublinear critical one Lum" (lu) oc Lu a occurs, as can be 
inferred from the inset of Fig.EIb). If T is far below T c , 
the knee disappears, and the spectrum is dominated by a 
single peak reflecting the AOP (see the next subsection). 
This holds for T = 0.8, 0.6, and 0.5. It is interesting to 
notice that there is only a weak temperature dependence 
in the position log w lu sa 1.5 of the high-frequency peak 
from high-T liquids down to deep-in-glass states. This 
feature will also be studied in the next subsection. 

Figurc[S]shows the reactive part of the normalized lon- 
gitudinal moduli A(uj) = 1 — um'(u>). Since m'(u>) and 
m"(ui) are related via the Kramers-Kronig relation [4J, 
features for A(u>) can be understood from those for the 
susceptibility spectra just mentioned. In particular, a 
peak in urn" {to) implies a strong increase of A(w) near 
that peak frequency. This way, one understands an al- 
most single-step increase of A(co) at log 10 cj w 1.5 for 
T 3> T c and T <C T c , reflecting the presence of only 
the high-frequency microscopic peak in com' 1 '(to) for those 
temperatures (c/. Fig-EJ. For T > T c , on the other hand, 
the presence of two peaks in Lum"(Lu) explains the two- 
step increase of A(lu) shown in Fig. E{a). Because of the 
strong T dependence of the a-peak frequency, the posi- 
tion where the low-w variation of A(lu) occurs depends 
strongly on the temperature. In addition, the sublinear 
susceptibility variation u)m"(u>) oc Lu a , present both in 
liquid and glass states near T c , also leads to the enhance- 
ment of A(lu). Indeed, it follows from Eq. I|26|) that 

A(lu) = 1 + C c + L>cos(7ra/2)r(l - a)(cut ) a , (29) 

holds near the (3 minimum for liquids or near the knee 
position for glasses. The curve based on this formula 
(dashed line) is compared with the memory kernel at the 
critical point (dotted line) in the inset of Fig.||jfb). 

For later convenience, let us define the a and (3 regimes 
for A(w). We start from liquid states. Since the memory 
kernel m(t) decays from the plateau C c to zero in the 
a regime, corresponding a regime for A(u>) shall be de- 
fined as the frequency interval < lu < lu^ in which lu^ 
satisfies A(w+) = 1 + C c . Thus, the a regime for A(lu) 
is below the horizontal bar marking 1 + C° in Fig. EJa). 
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FIG. 6: Reactive part of the normalized longitudinal moduli 
A(w) = 1 — wm'(w) for the memory kernels m(t) shown in 
Fig-EJfor liquid states (a) and for glass states (b). The insets 
exhibit the double logarithmic representation of A(u) for T 
close to T c . The meaning of the labels x = 2 and x = 3 is the 
same as in Fig. 0] Dotted line with label c denotes A(uS) at 
T c . Horizontal lines mark the plateau height 1 + C C (see text). 
Dashed lines in the upper inset show the MCT asymptotic- 
law results based on Eq. 130H for T — 1.8 and x — 3 (e < 0), 
while dashed line in the lower inset exhibits the curve based 
on Eq. 1291 . The open {ui^ ) and filled (uJp) diamonds refer 
to the frequencies uji defined in the text. 

In the p regime, the kernel m(t) is well described by the 
MCT asymptotic formula JUJ) (cf. Fig. 0. The corre- 
sponding [3 regime shall be defined as the one where A(w) 
can be well described by the asymptotic formula 

A(w) = l + C c -DwG'(w), (30) 

which follows from Eq. I|25|l . The dashed lines in Fig.EJa) 
show curves based on this formula. Let us introduce co^ 
as the frequencies where A(oj) (solid lines) differ from 
the MCT asymptotic formula (JSOJl (dashed lines) by 2%. 
These frequencies are marked by open (wj) and filled 
(k/g) diamonds in Fig. Ela), and the frequency interval 
Up < lu < u>p shall be considered as the (3 regime for 
A(lu). As mentioned in Sec. Ill CI there is an overlap 
between the a and (3 regimes. For glass states, the (3 
regime for A(lo) shall be defined as the frequency inter- 
val lu < uit where cut marks the point at which A(lu) 
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FIG. 7: (a) Solid line denotes the memory-kernel spectrum 
rh"(ui) for p = 1.093 and T — 0.5, and dashed and dotted 
lines respectively refer to its one-mode and two-mode contri- 
butions, m^ 1 - 1 (uj) and m' 2 ' (uj). Dash-dotted line shows the 
SGA spectrum rhg GA (uj) based on Eqs. I|31[l . (b) The same 
as in (a), but for T = 1.47. (c) Comparison of the memory- 
kernel spectrum m(uj) (solid line) with rh'^(uj) (dashed line) 
for the transversal current correlator and rh"(uj) (dotted line) 
for the velocity correlator for p — 1.093 and T = 0.5. 



differs from the MCT asymptotic formula by 2%. 
The frequency cu^ at the critical point is marked by the 
filled diamond in the inset of Fig. EJb) . 



C. Vibrational excitations 

In contrast to the structural-relaxation processes, there 
is only a weak temperature dependence in the time scale 
for the short-time dynamics in m(t) occurring at log 10 1 < 
— 1 (cf. Fig. 0J|. Correspondingly, the high-frequency 
peaks in the susceptibility spectra ujm"(uj) are located 
nearly at the same frequencies log 10 wrjI.5 from high-T 
liquid down to deep-in-glass state (cf. Fig.EJ), leading to 
a common increase of A(u>) around this frequency regime 
(cf. Fig.EI) • In this section, nature of such high-frequency 
microscopic process shall be discussed. 

We start from temperatures far below T c , T <C T c , 
called stiff-glass states, which are characterized by the 
Debye- Waller factors close to unity as exemplified by the 
curve for T = 0.5 shown in Fig.^ One understands from 
Eqs. JHU H9dfl. HUcjl . and (|10d|) that the renormalized 
coupling coefficients become very small for the stiff-glass 
states, and they decrease towards zero in the limit rj — 
1 — f q — > 0. Furthermore, the two-mode contributions 
to the kernels get suppressed relative to the one-mode 



contributions, in particular (t)/rh^ (t) — 0(rf). Fig- 
ure[7[a) shows the hatted spectrum rh"(ui) for T = 0.5 
(solid line) along with its one-mode rhS 1 ' (uj) (dashed 
line) and two-mode rh^ (uj) (dotted line) contributions 
(cf. Eq. JTOaJ). It is seen that the one-mode contribu- 
tion rrS 1 ^ (u>) indeed provides the dominant contribution 
to the total memory-kernel spectrum m"(u>). The role 
played by the two contributions rhS 1 ^ and rh^ is utterly 
different. The former describes the AOP, while the latter 
provides a background spectrum. 

Further properties of the AOP have been discussed 
in Ref. [24| based on a stiff-glass approximation (SGA), 
which is obtained by dropping the two-mode contribu- 
tions to the memory kernels. (See also Ref. [3(| on the 
SGA.) In particular, it has been found that harmonic os- 
cillations of the particles within their cages are a good 
description of the relevant modes, and the density corre- 
lator as well as the memory kernel in the stiff-glass states 
can be described as a superposition of independent har- 
monic oscillator spectra where the distribution of oscilla- 
tor frequencies is caused by the distribution of sizes and 
shapes of the cages. According to the cited MCT result 
for the AOP, one obtains m(ui) ~ ™sga(w) with 



wmsGA(w) = wi [ x(w) - 1 ] 



(31a) 



Here x(u>) is given by a superposition of undamped 
harmonic-oscillator spectra 



d£p(0xd"), 



(31b) 



m = VW-O (2^0, (31c) 

Xz(uj) = -SI 2 / [uj 2 - + iu)v\. (31d) 

The weight distribution for the oscillators with fre- 
quency extends from the low-frequency threshold 
f2_ = uj-VL to the high-frequency one S1+ = w+fi, where 
u>± = 1 ± y / tTii with an integrated coupling coefficient 

wi = ^Z fe (cf. Eq. l|lUcfl ). & denotes some aver- 
age value of £l q defined in Eq. Ij8b|l over the q range 
1 ^ <7d H3- Here, qu = (G^p) 1 / 3 denotes the Debye 
wave number, and qu = 4.02 for the density p — 1.093. 
Since the q range around the first peak position g max 
of the static structure factor S q is known to constitute 
the dominant contribution to the kernel 28], f2 in the 
present study is taken as an average of the first maxi- 
mum and first minimum of fl q , which occurs at q ~ qr> 
and q max , respectively (cf. Fig. I13fl . A friction term v is 
necessary when one wants to incorporate the neglected 
two- mode contributions in a pcrturbative manner |24| . 
but one may read the formulas with v — unless stated 
otherwise. Figure (7fa) demonstrates that Eqs. (f5T]) with 
v = (dash-dotted line) describe qualitative features of 
the AOP of m"(u>) for T = 0.50. 

By increasing T, one observes enhanced intensity in 
rh"(uj) at low frequencies as shown for T = 1.47 in 
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Fig. 0b) . This is a precursor phenomenon of the glass 
melting at the transition point. Near T c , such enhance- 
ment for decreasing frequencies can be described by a 
critical power law-w " 1 (cf. Eq. J2HI)- By substituting 
the critical decay of the density correlators (f> q (t) ~ t~ a 
into Eq. I|10b|) , one understands that it is the linear term 
rhP-'ft) that provides the critical decay of m(t), whereas 
the quadratic term ttj,^ 2 -* (t) serves only as a correction 
0(t~ 2a ). This explains why in Fig. Gfb) the strong in- 
crease of the intensity is dominated by (u>) (dashed 
line). Thus, the modes building the AOP are buried 
under the tail of the central peak, and they show up 
only as a shoulder. Nevertheless, it is remarkable that 
the SGA spectrum rhg GA (u>) (dash-dotted line) succeeds 
to extract the buried AOP portion as demonstrated in 
Fig. 0b). The shift of the AOP to lower frequencies upon 
increase of the temperature reflects the widened size of 
the cages so that particles "bounce" in their cages with 
smaller average frequencies. Indeed, the particles are lo- 
calized in such tight cages at T = 0.5 that the square 
root Sr — (Sr 2 (t — » oo)) of the long-time limit of the 
mean-squared displacement [3p| is only 5.0% of the par- 
ticle's LJ diameter, whereas it is increased to Sr — 0.139 
at T = 1.47. The shift is accompanied by an increase of 
the inelastic spectrum including the AOP, reflecting the 
decrease of the elastic contribution itC5(uj). 

The analysis presented in Fig. 0(b) also implies a pos- 
sibility to extract AOP portion in the spectrum on the 
basis of Eqs. I|31[) . So, let us examine to what extent 
the high-frequency microscopic peaks in the susceptibil- 
ity spectra Lom"(u>) at higher temperatures have such 
"harmonic" character. This is done in Fig. [S] where rep- 
resentative susceptibility spectra taken from Fig. are 
compared with the SGA susceptibility spectra 

wmg GA (w) = (l + C)«;ix"(w) J (32) 

which follows from Eqs. (|%c|) and H31a(l . To focus on the 
high-frequency regime, the comparison is done in Fig. 
with the linear-w axis. In evaluating the SGA susceptibil- 
ity spectra for liquids, nonergodicity parameters f q and 
C in the equations and quantities involved are replaced 
by the ones at the critical point. 

As discussed in connection with Figs. [3(a) and 0b) , 
the SGA spectra describe the AOP portion of the 
memory- kernel spectra for T = 0.5 and T = 1.47 fairly 
well, and this is reflected in the susceptibility spectra 
shown in the two bottom panels of Fig. Somewhat 
large deviation seen at w > 40 for T = 1.47 is due to 
the neglected two-mode contribution, which is more en- 
hanced than the deviation in Fig. 0b) due to the presence 
of the factor to in ujm"(uj). On the other hand, because 
of this factor, the deviation at small frequencies seen in 
Fig. 0b) is suppressed for turn" {to). 

One understands from the other panels of Fig. that 
the harmonic approximation based on Eq. I|32l) reason- 
ably accounts for overall features - peak position and 
strength - of the microscopic peaks in ujm"{u>) not only 
near T c but even at high temperatures such as T = 10. 
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FIG. 8: Susceptibility spectra um"(u) (solid lines) at indi- 
cated temperatures taken from Fig. but plotted here on 
the linear-tj axis. Dashed lines denote the SGA susceptibility 
spectra uttig GA (ij) based on Eq. 1321 . 



The large deviations for low frequencies at T = 4.0 
and 10 simply reflect the fact that at such high T the 
structural-relaxation contributions also enter into the 
high-frequency regime (cf. Fig. 0). Thus, considerable 
portion of the short-time or high-frequency microscopic 
process, from high-T liquids down to deep-in-glass states, 
can be described as a superposition of harmonic vibra- 
tional dynamics. This explains the weak T dependence 
of the microscopic process shown in Figs. 14161 

The persistence of the vibrational dynamics inside the 
cage in the liquid state is not surprising. To see this, we 
show in Fig. the velocity correlator ^(t) for represen- 
tative liquid and glass states. It is well known that the 
cage effect in dense liquids manifests itself by oscillatory 
variations with a decay of ^>(t) to negative values [j], 
and Fig. EJa) demonstrates this phenomenon. That such 
oscillatory dynamics exhibited by reflects the AOP 
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FIG. 9: Normalized velocity correlator at p = 1.093 

for temperatures referring to liquid states (a) and glass states 
(b). 



can be understood from the observations (i) in the stiff- 
glass state T = 0.5 the spectrum m"(w) of the relaxation 
kernel for ^(t) exhibits an AOP whose spectral features - 
threshold frequencies and maximum position and height 
- nearly coincide with those for the AOP of rh"(uS) as 
demonstrated in Fig.[7fc), and (ii) at higher T including 
liquid states there is a correlation between the peak posi- 
tion in the SGA susceptibility spectra tvmg GA (u>) shown 
in Fig. [H] and the time at which the velocity correlator 
'F(i) reaches the first minimum. 

There is another interesting feature caused by the 
dominant harmonic nature of the dynamics in stiff-glass 
states. From A(w) for T — 0.5 shown in Fig. EJb), one 
observes the decrease of A(w) from A(ui = 0) = 1 + C 
at log 10 w ~ 1. One can show that such behavior is also 
well described by the harmonic approximation: it fol- 
lows from Eqs. (|Sc|) and l|31a|l that the reactive part of 
the modulus within the SGA is given by 

AsgaM = (1 + C) { 1 - Wl [x'M - 1 ] }. (33) 

Figure EH compares A(w) for T = 0.5 with the curve 
based on Eq. i|33|) on the linear-w axis. One understands 
that decrease of A(w) from A(w = 0) = l + C observable 
at lo w 10 is due to the harmonic modes near the thresh- 
old frequency f2_. Similarly the bump at w ss 60 reflects 
the harmonic modes near Q + . 
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FIG. 10: Reactive part of the longitudinal modulus A(ui) 
(solid line) for p = 1.093 and T = 0.5 taken from Fig. El but 
plotted here on the linear-o; axis. Dashed line exhibits the 
SGA modulus Asga(w) based on Eq. 1331 . Horizontal dotted 
line denotes A(u> = 0) = 1 + C. Vertical arrows mark the 
threshold frequencies f2_ and fi+ defined in the text. 



D. Structural-relaxation and vibrational-dynamics 
contributions to sound-velocity dispersion 

Here we provide an explanation of the anomalous 
sound- velocity dispersion shown in Figs. |21 and [21 based 
on Eq. (|24bjl and on the features of A(uj) explored in the 
previous two subsections. 

At T = 10, the whole relaxation of the memory kernel 
m(t) occurs on the short time scale log 10 t < — 1 (Fig.0J. 
This yields a susceptibility spectrum uira"(ui) consisting 
only of the microscopic peak located at log 10 w « 1.5 
(Fig. 03, and the entire increase of A(w) occurs in the 
frequency regime log 10 tj > (Fig. EJ. It is therefore 
sufficient to use the linear-w axis to describe the whole 
lo variation in A(lo), and this explains via Eq. I|24b|) the 
observation in Fig.|2a) that the entire positive dispersion 
of the sound velocity v q starting from vq up to close to 
Voo occurs within the linear-g axis. 

By decreasing the temperature to T = 4.0, the relax- 
ation of m(t) becomes slower and stretched due to the 
development of the cage effect (Fig. 0J, leading to the 
evolution of low- frequency contributions in torn" (oj) and 
A(w) (Figs. |U and EJ. Translated to v q via Eq. J24bj, 
this means that the approach of v q towards the hydro- 
dynamic value vq becomes delayed compared to the one 
at T = 10. This explains why v q for small wave numbers 
shown in Fig. Efb) is still located above vq. 

For temperatures close to but above T c , the mem- 
ory kernel m(t) exhibits the glassy structural relaxations 
(Fig. 2}, leading to low-u; variations in ujm"(uj) and A(w) 
which can be adequately described only with the log 10 ijJ 
axis (Figs. [5] and 0). Let us analyze the sound- velocity 
increase for T = 1.8 in detail. For this purpose, we replot 
the inset of Fig. EJc) in Fig. IllT ah but with the log 10 Q 
axis. This is necessary to fully appreciate the small-g 
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FIG. 11: (a) Apparent sound velocity v q = u\ max /g (circles) 
for p — 1.093 and T = 1.8, taken from the inset of Fig. ffic). 
but plotted here with the log 10 q axis. Long-dashed line is 
from the solution to Eq. 12311 . while solid line is based on 
Eq. (|24bpi . Horizontal lines denote vo (dashed line), Vao (dash- 
dotted line), and (dotted line). The contributions Av a , 
Av/3, and Av^ to the sound- velocity increase defined in the 
text are indicated by vertical arrows. The open and filled 
diamonds respectively mark and defined in the text, 
(b) The same as in (a), but here the sound velocity refers 
to T — T c , taken from the inset of Fig. |^a), and horizontal 
dashed line denotes vq = uoVl + C. 



variation of the sound velocity v q , which on the basis of 
Eq. (I24bl) originates from the glassy low-tj variations in 
A(w). Notice that the circles from the full MCT solutions 
for (fir (lu) can be obtained only up to the minimum value 
of the discretized wave-number grids. On the other hand, 
there is no such restriction in the solution (solid line) to 
the lowest-order equation (|24bjl of the generalized hydro- 
dynamic description. The difference between the circles 
and the solid line - in the q range where the circles are 
available - is due to the approximations involved, but 
there should be no problem to use the approximate solu- 
tion for understanding the essence of the anomalies. 

Let us quantify the contributions to the sound velocity 
increase due to the a process (to be denoted as Av a ), 
due to the (3 process (Avp), and due to the microscopic 
process (Av^) in the following way based on the a and 
(3 regimes for A(w) defined at the end of Sec. IIII Bl As 



mentioned there, A(u) varies from A(ui = 0) = 1 to 
A(w+) = 1 + C c in the a regime. One therefore obtains 
from Eq. (|24E|l : Av a = u£, - v with = v y/l + C c . 
The notation comes from the fact that v^ is the high- 
frequency limit of the sound velocity in the a regime. 
The horizontal dotted line in Fig. Illf a) marks v^, and 
corresponding linear dispersion law qv^ and the veloc- 
ity have been included with dotted lines in Fig. |2Jc). 
The /3-relaxation contribution Avp can be obtained in 
terms of the frequencies ui^ also introduced at the end 



of Sec. 1111 Bl one obtains from Eq. (|24bfl A 



V 



u 



with v^ = vq^Ja(lj^). These velocities are marked by 

open (vp) and filled (v~p) diamonds in Fig. ITlT a). The 
rest of the sound-velocity increase towards Voa is due to 
the microscopic process. As mentioned in the previous 
subsection, considerable portion of the microscopic pro- 
cess is due to the harmonic vibrational dynamics. The 
contributions Av a , Avp, and Av^ are indicated by verti- 
cal arrows in Fig. lllf a). Since there is an overlap between 
the a and (3 regimes (cf. Sec. Ill C|l . there is a correspond- 
ing overlap between Av a and Avp. 

The full sound-velocity variation for T — 1.8 can be 
understood in this way in terms of the contributions 
from a- and /3-relaxation processes and from the micro- 
scopic vibrational dynamics. Notice that the structural- 
relaxation contributions can be adequately described 
only with the logarithmic q axis. This explains why in 
Fig- EIc) the approach of v q towards the hydrodynamic 
value vq could not be resolved, but only the approach 
towards v^ can be observed. Thus, the sound-velocity 
variation shown in Fig. |2Ic) mostly reflects the micro- 
scopic contribution Ai^. 

In glass states, the a process is arrested, and the hy- 
drodynamic sound velocity is altered to vq = t>o\/l + C. 
There are at most f3 and \x contributions to the sound- 
velocity increase from this hydrodynamic value. At and 



near T C: the (3 contribution, Avp 
neglected, where v~£ 



Vq, cannot be 



vq^J A(u}p) with ujp for glasses 

defined at the end of Sec. IIII Bl This feature can be 
understood from Fig. Illf b). a redrawing of the inset of 
Fig-EHa) with the logarithmic q axis, where the (3 contri- 
bution Avp and the microscopic contribution Av^ (de- 
fined as the rest of the sound-velocity increase towards 
v<x>) are marked by vertical arrows. The /3-relaxation 
contribution Avp cannot be resolved with the linear-g 
axis, and this explains why in Fig. |3f a) the approach of 
v q towards vq with decreasing q is retarded. 

For temperatures further away from T c , the (3- 
relaxation contribution becomes smaller, and the sound- 
velocity increase from vq is entirely due to the micro- 
scopic process. This explains why in Fig. OJb) the ap- 
proach of v q towards vo can be observed for T — 1.47 
even with the linear-g axis. By further decreasing the 
temperature to T = 0.5, a new feature shows up reflect- 
ing the dominant harmonic nature in the stiff-glass-state 
dynamics. As discussed in connection with Fig. 1101 there 



15 



appears a frequency interval where A(o>) becomes smaller 
than A(u) = 0). According to Eq. (|24b(l . this gives rise 
to a "negative" dispersion, and explains the result for 
T = 0.5 shown in Fig. Etc). 



IV. LOW-FREQUENCY EXCITATIONS IN 
LONGITUDINAL AND TRANSVERSAL 
CURRENT SPECTRA 

The evolution of the AOP is related to the arrest of 
density fluctuations caused by the cage effect, and the 
AOP is the result of a mapping of the cage distribution 
on the frequency axis (c/. Sec. lIIICT) . Since the arrest of 
density fluctuations is driven by the ones with a wave 
number q near the structure- factor-peak position |28j . 
this interpretation of the AOP suggests that it appears 
in spectra of all probing variables that couple to density 
fluctuations of short wavelengths. But different probing 
variables will weight the oscillating complexes differently, 
and therefore the shape and the peak position of the AOP 
will depend somewhat on the probe. Two such exam- 
ples have been analyzed in Ref. 2A\: the tagged-particle 
density correlators and the velocity correlator In 
particular, it has been shown there how the AOP in the 
memory-kernel spectrum m"(w) for (cf. Fig. Etc)) 
accounts for the excess intensity in the density of states 
2^"(u))/tt with respect to the Debye model. 

In this section, we investigate how the AOP manifests 
itself in spectral features of the longitudinal and transver- 
sal current spectra. This issue is relevant since in re- 
cent computational studies H, 0] two excitations have 
been observed in the longitudinal current spectra, but 
the low-frequency boson-peak-like excitations have been 
assigned to the transversal mode. Since there is no direct 
cross correlation between the longitudinal and transver- 
sal current dynamics, the appearance of the transversal 
mode in the longitudinal spectra is interpreted as being 
due to "mixing" phenomena [lsf. This is in contrast to 
the above interpretation of the AOP, according to which 
possible low-frequency excitations in both of the longitu- 
dinal and transversal current spectra should be the man- 
ifestation of the AOP. Indeed, the AOP is present also 
in the memory kernel for the transversal current corre- 
lators, and its spectral features are quite similar to the 
one for the longitudinal current correlators as can be in- 
ferred from Fig. Etc). To corroborate our statement, we 
will also study the density dependence of spectral fea- 
tures, because such study was claimed to support the 
interpretation in terms of the mixing of the longitudinal 
and transversal modes (c/. Sec.[|J. 

A. Manifestation of AOP in spectral features 

We present with solid lines in Fig. ^] typical density 
fluctuation spectrum <p'^{uj), longitudinal current spec- 
trum <fi^ {ui), and transversal current spectrum (f^ (u) 
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FIG. 12: (a) Density fluctuation spectrum <f>'^(iv), (b) longi- 
tudinal current spectrum <f>]j (u>), and (c) transversal current 

spectrum (f^ (w) at wave number q = 3.7 for p = 1.093 and 
T = 0.5. Solid lines are based on the full MCT equations of 
motion described in Sec. [H] Dash-dotted lines show curves 
within the generalized hydrodynamic description, Eqs. 134H . 
Dashed lines denote the HFS part based on Eqs. 13511 . whereas 
dotted lines refer to the AOP portion from Eqs. 113611 . Curves 
for 4>q(uj) are obtained from those for 4>\ ipj) via Eq. 1131 . 

at deep in the glass state, T = 0.5 <C T c . These results 
are based on the full MCT equations of motion described 
in Sec. [H] The wave number q — 3.7 chosen for the 
demonstration is about half of the structure-factor-peak 
position (c/. Fig. The fluctuation spectrum <j>g(u>) 
shown in Fig. I12f a~) exhibits two peaks for ui > 10, and 
illustrates a hybridization of the high-frequency sound 
with the modes building the AOP [24J : the narrow peak 
located at higher frequency is due to the high-frequency 
sound propagation, and a broad lower-frequency excita- 
tion reflects the AOP studied in Sec. IIII CI On the other 
hand, there is apparently only one peak in the longitudi- 
nal current spectrum (p^ (oj) shown in Fig. I12r bh It is 
due to the high-frequency sound, and the peak position 
w Lmax j g near iy the same as the one for the fluctuation 
spectrum. It is not clear how the AOP manifests itself in 
the spectral shape of <fi^ (u>) since, compared to <j)q(u), 
the low frequency part is suppressed due to the presence 
of the factor ui 2 (cf. Eq. (|13|) h The transversal cur- 
rent spectrum ^ (uj) shown in Fig. Et c ) a l so exhibits 
only one peak. Its peak frequency wj max is located at a 
lower frequency than ui^ max reflecting a smaller transver- 
sal sound velocity ojJ ma,K /q compared to the longitudinal 
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FIG. 14: (a) Solid line denotes the memory-kernel spectrum 
m(u>) for p = 1.093 and T = 0.5, while dash-dotted line refers 
to the SGA spectrum based on Eqs. 1311 . Arrows mark the 
threshold frequencies £l± defined below Eqs. I31H . Dotted line 
denotes the frequency dependence of the function <fif AOP (uj) 
given in Eq. I36at in arbitrary units, and its peak frequency 
waop is marked by an arrow. Long dashed line denotes the 
corresponding result based on Eq. 137aH . (b) The same as in 
(a), but for p = 0.95 H3]. 



FIG. 13: Bare dispersion relation Q q for the longitudinal 
sound (solid line) and Cl^ for the transversal sound (dashed 
line) at T = 0.5 for p = 1.093 (a) and p = 0.95 (b). Hori- 
zontal dotted lines denote the threshold frequencies £l± (cf. 
Fig. 1141 , and arrows mark the Debye wave number qd and 
the structure- factor-peak position q max . 



frequencies are marked in Fig. I14f a). The hybridization 
occurs in the range fi_ < fl q < f2+ indicated in Fig. H3f a) 
where the two oscillations overlap. The generalized hy- 
drodynamic description for the transversal current can 
be introduced similarly on the basis of Eq. (|18fl : 



one. Again, it is not clear whether there is a contribution 
to the spectral shape of (f^ (uj) from the AOP. It is the 
purpose of this subsection to clarify this point. 

To this end, we will derive approximate formulas spe- 
cialized to describe each of the high-frequency-sound 
(HFS) and low-frequency-AOP portions of the spectra. 
The starting point is the following expression under the 
generalized hydrodynamic description (cf. Eq. (1221 and 
citation [23|), obtained from Eq. Ijl4|l with the hatted 
memory kernel rh q (ui) approximated by its q — * limit, 
m(u>) = m q=0 (uj): 



uj 2 Cl 2 q m"(uj) 



[uj 2 - n 2 q A(uj)} 2 + [ujn 2 q m"(uj)f 



(34a) 



Here A(w) = 1 — ujrh'(uj). Corresponding description for 
the fluctuation spectrum (j) q (uj) is derived via Eq. (fLl)) . 
The formula (|34a[) describes hybridization of two oscil- 



lations: one is a bare sound with dispersion whose 
wave-number dependence is shown in Fig. ll3f a). and the 
other is represented by the AOP spectrum in rh"(u>). The 
AOP region in the spectrum m"(u>) shall be defined as 
Sl_ < uj < tt + in terms of the threshold frequencies £l± 
introduced via the SGA spectrum in Sec. IIII CI These 



uj 2 (nJ) 2 m^(Lu) 



[UJ 2 - (^) 2 A T (uj)} 2 + [ W (fiT)2^4( W )]S 



(34b) 

mj_ (w) and At(w) = 1 — ujrh' T (uj). 
rh'^(uj) as shown in Fig. [7f c), one un- 



Here wit(w) 
Since m"(uj) 

derstands that the hybridization of the bare transversal 
sound of dispersion Cl q with the AOP occurs also in the 



range £7_ < fl q < Sl + indicated in Fig. I13f a). Fig- 



ure ^] demonstrates that the generalized hydrodynamic 
description (dash-dotted lines) reproduce main features 
of all spectra fairly well. In particular, the subtle hy- 
bridization of the high-frequency sound and the AOP for 
4>q(uj) is treated semiquantitatively correctly, justifying 
the use of Eqs. I|34l) as starting equations. 

The HFS portion of (f>^ (uj) shall be approximated by a 
Lorenzian (multiplied by uj 2 ) with the width determined 
at the peak frequency w^ max : 



,L-HFS" 



L0 2 n 2 m"(uj^ max ) 



( w L ma x)2]2 



[ujVl 2 q m" (u)\ ma *)] 2 



(35a) 

The dashed curve in Fig. I12f b) shows the result based 
on this formula, and the corresponding HFS curve for 
the fluctuation spectrum oc 0^~ HFS (u))/uj 2 is included in 
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Fig. I12f a') . One understands that the formula l|35a|) well 
describes the HFS portion of these spectra. Thus, the 
rest of the spectral intensity in the low-frequency regime 
should be due to the AOP. 

Corresponding expression for the HFS part of the 
transversal current spectrum reads 



^T-HFS" 



1 ^ 2 m^(wJ max 



( w Tmax)2]2 



[u/(QT)2^//( w Tm ! «)]2 

(35b) 

Figure Et c ) shows that this formula (dashed line) fairly 
catches the main peak of the spectrum, but the predicted 
width is considerably narrower than the one of the full 
MCT solution (solid line). This feature can be under- 
stood in terms of the AOP contribution as follows. We 
first notice that wj max « 39 for q = 3.7 is located in 
the central region of the AOP (cf. Fig. I14f a1'). Let us 
introduces frequencies <D± in this central regime so that 
u)_ < cjj max < uj + . Since there holds m^(cj) ps rh"(Lo) as 
shown in Fig. Etc), the Kramers-Kronig relation implies 
At(w) « A(w). One therefore understands from Fig. I1UI 
that At(w) in the central regime of the AOP increases 
with u, i.e., A T (w_) < A T (wJ max ) < A T (w+). Since 
in lowest order wj max is the solution to the equation 
( w Tmax)2 = (QT) 2 A T (oj^ max ) (cf. Eq. J23J), one has 

Tmax\2 



(Qj) 2 A T (£-) _< (a,, 
with 



< (ttJ) 2 A T (£ 4 



Combined 



_ < (u>g max ) 2 < this leads to an inequality 



\Q 2 (> <t 



I . , 1' max 



J ± " V^g )] 2 ^ L^± , 

tice from Fig. Etc) that there is only a weak oj dependence 
in m'^(u)) in the central regime of the AOP. All this to- 
gether, one finds on the basis of Eqs. (j34b|) and (|35b(l : 

(fil (<^±) > < /'J~ HFS (<^±)- This way, one understands 
that the residual intensity in Fig. I12f c1 which cannot be 
described as the HFS portion (dashed line) is caused by 
the frequency dependence of At(lo) near wj max , which 
in turn is due to the AOP. From the mentioned reason- 
ing, it is clear that such deviation from the Lorenzian 
shape of the spectrum occurs whenever the sound reso- 
nance frequency is located in the central regime of the 
AOP. This explains why such deviation is small for the 
longitudinal current spectrum shown in Fig. H2f b). whose 
resonance frequency ui^ max ps 65 exceeds f2+ and is not 
located in the central regime of the AOP. Such larger 
resonance frequency max than expected from the bare 
dispersion Cl q for q = 3.7 shown in Fig. I13f a~) is due to 
the positive dispersion discussed in Sec. IIIII 

In view of the results presented so far, one understands 
that there are residual spectral intensities due to the 
AOP in the frequency regime to < ujq^ max . In the fol- 
lowing, we will derive approximate formulas extracting 
such AOP portion based on Eqs. I|34|l . For this purpose, 
let us consider the denominator on the right-hand side of 
Eq. (|3lajl . [w 2 - Cl 2 q A(w)] 2 + [u f2 2 m"(w)] 2 . The appear- 
ance of the HFS peak is mainly due to the decrease of the 
first term towards zero approaching the resonance posi- 
tion t^ max . Since we are interested in the off-resonant 



(ftJ) 2 A T (w±)l 2 . We also no- 



regime io < w^ max , this term will be approximated by 



fi*. Corre- 



its the small-w limit, [to 2 — Cl 2 A(uj)] 2 
spondingly, only the leading-order contribution for small 
uj shall be retained in the second term, [u> ft 2 m"(ui) ] 2 « 
[ v & 2 rh"(uj = 0) ] 2 . This results in the following expres- 
sion for possible AOP portion in (fi^ (uS): 



io 2 m"(io) 



Q 2 1 



[wm"(uj = 0)} 2 



(36a) 



According to the derived formula, the q and to depen- 
dences are factorized. This implies that the peak position 
of the AOP portion - when it is visible - is q independent. 
The frequency dependence of the formula (|36ajl is shown 
as dotted line in Fig. 114( a). and its peak frequency, to 
be denoted as waop? is marked by an arrow. Further- 
more, one infers from the q dependence of Cl q shown in 
Fig. ll3l that the intensity of the AOP portion is predicted 
to increase upon increase of the wave number for q < grj, 
whereas this trend becomes reversed for q }t qr>. Cor- 
responding formula for the transversal current spectrum 
can be derived starting from Eq. I(34b|l : 



6 T-AOP" 



' 2 774(w) 



(fiT)2 1 + [cum^LO = 0)] 5 



(36b) 



Since rh^(w) « rh"(u) (cf. Fig. Etc)), Eqs. © imply 
that the peak positions of the AOP portions of both lon- 
gitudinal and transversal current spectra - when they are 
visible - are located at the same frequency waop- Fur- 
thermore, a stronger intensity of the AOP portion is pre- 
dicted for the transversal current spectra in the pseudo 
first Brillouin zone since there holds Ct q > ClJ as can be 
seen from Fig. I13f a) . Finally, we notice that it is reason- 
able to further approximate the above formulas as 



6 L-AOP" (a;) = 2 



uj 2 m^ 



& q 1 + [wm( 2 )"(w = 0)] 5 



6 t-aop" 



(^) 2 l + [com^"(io = 0)f 



(37a) 
(37b) 



since it is the one-mode contribution which provides 
the dominant contribution to the memory kernel and 
which describes the AOP, whereas the two-mode con- 
tribution provides a background spectrum as discussed 
in Sec. IIII CI Figure Et a ) demonstrates that indeed the 
difference between the curves based on Eqs. (|36H and l|37|l 
is small. Our motivation for introducing Eqs. I|37|l will 
become clear in the next subsection. 

The dotted lines in Fig. El are based on Eqs. iJSSJ. For 
the density fluctuation spectrum (Fig. I12f a)). it is seen 
that the dotted line, obtained from Eq. (|36a|l via the rela- 
tion 113(1 . extracts the AOP portion quite reasonably. It 
is also seen from Fig. ll2f b') that Eq. I|36a|) provides a rea- 
sonable account of the residual low-frequency intensity in 
the longitudinal current spectrum located at w « waop- 
Notice that the AOP appears differently in <j) q (oS) and 
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FIG. 15: Longitudinal current spectra 0, (w) (solid lines) at 
p — 1.093 and T = 0.5 for indicated wave numbers. Dashed 
lines denote the HFS part based on Eq. I|35a|l . whereas dotted 
lines refer to the AOP portion obtained from Eq. I36at . 



p = 1 .093 




FIG. 16: Transversal current spectra </>J (lo) (solid lines) at 
p — 1.093 and T = 0.5 for indicated wave numbers. Dashed 
lines denote the HFS part based on Eq. d35bll . whereas dotted 
lines refer to the AOP portion obtained from Eq. I36bl . 



cj)^ (lu) reflecting the difference of the factor lo 2 . To un- 
derstand the result shown in Fig. ll2f c) for the transversal 
current spectrum, we notice that, by construction, the 
applicability of the formula l|36b(l is limited to the small- 
lo regime considerably lower than the sound resonance 
position wj max . Therefore, the dotted line in Fig. IP2T c) 
in the region lo > wj max has no physical meaning. Nev- 
ertheless, it is seen that the dotted line accounts for the 
residual intensity in the low-frequency regime which can- 
not be described as the HFS portion (dashed line). Thus, 
Eq. (|36b|) describes the low-frequency part of the devia- 
tion from the Lorenzian shape explained above. 

The variation of the longitudinal and transversal cur- 
rent spectra (solid lines) with changes of wave number q 
is presented in Figs. Upl and 1 161 along with their decom- 



position into the HFS part (dashed lines) and possible 
AOP portion (dotted lines). The dispersion relation and 
the full width at half maximum (FWHM) as a function 
of q are summarized in Fig. 1171 We start from the dis- 
cussion on the longitudinal current spectra. When q is 
small so that Wg max < 17-, there is practically no ef- 
fect from the AOP on the spectral shape. (However, 
the AOP affects the resonance frequency w^ max near fi_ 
as discussed in Sec. IIII Dl ) Therefore, the spectrum for 
q = 0.5 shown in Fig. ^|is dominated by the HFS por- 
tion (dashed line), and there is no physical meaning in 
the dotted line for this q value. Thus, the spectrum for 
q = 0.5 is very close to Lorenzian (multiplied by to 2 ), and 
its width agrees with the one from the hydrodynamic pre- 
diction q 2 VQm"(uj = 0) as demonstrated in Fig. llTf b). As 
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FIG. 17: (a) Peak positions c^ max 



of 



and w, 



T max 



of 



g (w) (open circles) 
.i (a;) (filled circles) as a function of q for 
p = 1.093 and T — 0.5. Solid and dashed lines denote the bare 
dispersions Cl q and Cl£ , respectively. Dotted lines exhibit the 
hydrodynamic dispersion laws qvo and qvj . Peak positions of 
the AOP portion in <p q (cu) are denoted by open diamonds. 
Horizontal arrows mark Q± and ojaop taken from Fig. I14f a'l . 
(b) Full width at half maximum (FWHM) for <j>\" (w) (circles) 
as a function of q. Dashed line denotes the Lorenzian width 
Q.grh"(ujq max ) predicted by Eq. <l3"5al . Dotted line refers to 
the hydrodynamic width, q 2 VQifi"(uj = 0). Vertical arrow 
marks the wave number where c^ max = holds, (c) The 
same as in (b), but for <^ (w) for which the Lorenzian width 
(dashed line) is given by (f7^) 2 mT(wJ max ) and the hydrody- 
namic width (dotted line) by q 2 (vo) 2 rh'^(u> — 0). Vertical 
arrow marks the wave number where u> J max = Q.- holds. 



q increases so that max becomes larger than J7_ , effects 
from the AOP set in. These effects increase the spectral 
width, as exemplified for q = 1.3 and 2.1 in Fig. 1151 in two 
ways. First, the width becomes larger compared to the 
hydrodynamic width because tji"(uj^ max ) > m"(u> = 0) 
when the resonance frequency u^ max is located in the 
AOP regime, fi_ < ^ max < fl + (c/. Fig.H¥a)). This 



shows up as the deviation of the dashed line from the 
dotted line in Fig. I17f b). In addition, the frequency de- 
pendence of A(u>), which again is caused by the AOP, 
also enlarges the width as explained above in connec- 
tion with Fig. I12f c). This latter effect leads to a the 
non-Lorenzian shape of the spectrum, i.e., the deviation 
of the solid line from the dashed line for q = 1.3 and 2.1 
shown in Fig. 1151 and also explains the difference between 
corresponding circles and dashed line in Fig. I17f b) . The 
dotted lines in Fig. 1151 account for the low-frequency part 
of the non-Lorenzian spectra as discussed above. As q is 
increased further, the positive dispersion effect becomes 
important, and it pushes o^ max considerably above fre- 
quencies expected from the bare dispersion Cl q as shown 
in Fig. I17f a). Thereby, there opens a frequency window, 
where waop < ^\ max holds and the appearance of the 
AOP itself can be observed in the spectra. This fea- 
ture holds q > 2.9 as shown in Fig. ^| As mentioned 
above, there is no q dependence in waop> which results 
in flat dispersion curve for the low-frequency excitations 
as demonstrated with open diamonds in Fig. I17t a). 

Corresponding results for the transversal current spec- 
tra shown in Figs . IT61 and 1 1 71 can be explained similarly. 
But here, the positive dispersion of wj max from the hy- 
drodynamic law qiiQ with Vq 



Cj =0 , obtained from 



the leading-order expansion of Cl^ , is not so strong as in 
the longitudinal case. This is because the bare transver- 
sal sound dispersion f2j for q < qu is located below 
the maximum frequency u) w 40 of the the susceptibility 
spectrum ujrh"{u>) (cf. the bottom panel for T = 0.5 in 
Fig-El), where the positive dispersion effect is most signifi- 
cant as explained in Sec. lIIII Thus, there is no frequency 
regime where cjaop < wj max holds for the transversal 
current spectra, which can be inferred from Fig. I17t a). 
Therefore, only the width of </>J (ui) is affected by the 
AOP, which is summarized in Fig. I17f c). However, this 
feature is altered when the density is decreased, as we 
will see in the next subsection. 



B. Longitudinal and transversal current spectra at 
lower density 

In this subsection, effects of varying density on spectral 
features of the longitudinal and transversal current spec- 
tra shall be investigated. Specifically, a density p = 0.95 
lower than 1.093 studied so far will be considered, but 
with temperature T = 0.5 fixed. 

Upon decrease of the density, the bare sound dis- 
persions Cl q and and the memory-kernel spectrum 
rh"(uj) are shifted to lower frequencies as demonstrated 
in Figs. ll3f b) and !14f b). respectively. Notice the striking 
similarity between the change in m"(u>) upon decrease of 
the density at fixed T and the change upon increase of 
the temperature at fixed p found in Fig. Efb). Indeed, 
the critical temperature is decreased to T c w 0.644 for 
p = 0.95, and lowering the density at fixed T = 0.5 drives 
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FIG. 18: The same as in Fig. 1151 but for a lower density 
p = 0.95, and dotted lines referring to the AOP portion are 
based on Eq. I37at as discussed in text. 
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FIG. 19: The same as in Fig. 1161 but for a lower density 
p = 0.95, and dotted lines referring to the AOP portion are 
based on Eq. 13 7b0 as discussed in text. 



the system closer to the critical point. As discussed with 
Fig. |7Jb) , a strong central peak is formed when the sys- 
tem approaches the critical point, and the modes building 
the AOP get buried under its tail. But again, Fig. I14f b) 
demonstrates that the SGA kernel room (lo) extracts the 
buried AOP portion quite reasonably ^37] . Thus, in anal- 
ogy to what is found in Fig.[7Jb), the AOP shifts towards 
lower frequencies upon decrease of the density, accompa- 
nied by an increase of its spectral intensity. 

As a related problem, we found that the use of 
Eqs. (|36[) to extract the AOP portion from current spec- 
tra leads to unplausible results. This is because of the 
mentioned development of the central peak, leading to 
quite large value of m"(u> = 0). As will be demonstrated 
in the following (cf. Figs. I^Of b) and I2l)f c)). such large 
value of m" (lo — 0) caused by the critical decay near T c is 



not appropriate for handling the high-frequency dynam- 
ics. We will therefore use Eqs. (|37|) instead of Eqs. 136( 1. 
which essentially amounts to replacing m!' (lo = 0) with 
rh^ (to = 0). The curves based on Eqs. (|S7|) arc not 
much affected by the critical dynamics because (i) only 
corrections to the critical decay enter into (lo) as ex- 
plained in connection with Fig.JTJb), and (ii) the critical- 
decay contribution in rrv 1 ' (lo) is suppressed in Eq. I|37a|) 
due to the presence of the factor lo 2 in its denomina- 
tor. As demonstrated in Fig. I14f ah there is practically 
no difference between the use of Eqs. <|36|) and (|37|l for 
p = 1.093, and we expect that physics is not altered with 
the use of the latter for p = 0.95. 

Longitudinal and transversal current spectra for p = 
0.95 (solid lines) with changes of wave number q are pre- 
sented in Figs. 1181 and !19l along with their decomposition 
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at smaller q for p — 0.95, as can be inferred by comparing 
Figs. 1171 and 121)1 These features are also due to the shift 
of the AOP to lower frequencies upon decrease of the 
density. Furthermore, there is no "negative" dispersion 
(c/. Sec. IIII Dl for p = 0.95 since the memory kernel for 
T = 0.5 at this density is not dominated by the harmonic 
contributions (cf. Fig. ED )))- Otherwise, the variation 
of the longitudinal current spectra with changes of q for 
p = 0.95 is quite similar to the one found for p = 1.093, 
and the discussion shall not be repeated. 

The variation of the transversal current spectra with 
changes of q for p = 0.95 is also quite similar to the 
one found for p — 1.093, as can be understood by com- 
paring Figs. El an d El But here, the resonance fre- 
quency uij max of the transversal sound exhibits a strong 
positive dispersion as shown in Fig. I20f a). Again, this 
is caused by the shift of the AOP to lower frequencies. 
The positive dispersion of the transversal sound velocity 
due to the contribution from the AOP can be discussed 
quite similarly to what is presented in Sec. IIIII for the 
longitudinal sound. Thereby, there opens a frequency 
window also for the transversal current spectra, where 



waop < w, 



T max 



holds and the appearance of the AOP 



itself can be observed in the spectra. Such AOP portion 
can be observed for q > 3.7 in Fig. El yielding a flat 
dispersion curve for the low-frequency excitations of the 
transversal current spectra as demonstrated with filled 
diamonds in Fig. I20f a). This is in contrast to the result 
found for p = 1.093, where u>aop > wj max holds for the 
whole q range and the effects from the AOP showed up 
only in the width of the transversal current spectra. 



FIG. 20: The same as in Fig. 1171 but for a lower density 
p = 0.95. In (a), filled diamonds are added denoting the 
peak positions of the AOP portion in ci^ (lo). Dash-dotted 
lines are added in (b) and (c) denoting q 2 VQrhP^ (oj — 0) and 
q 2 (vQ) 2 m^ {lu = 0), respectively (see text for details). 



into the HFS part (dashed lines) and possible AOP por- 
tion (dotted lines). As mentioned above, the dotted lines 
in these figures are based on Eqs. I|37|l . The dispersion 
relation and the FWHM as a function of q are summa- 
rized in Fig. El It is seen from Figs. I2()f b) and I20f c) 
that indeed the hydrodynamic prediction for the width 
with rh"(u> = 0) does not work for p = 0.95; instead, the 
corresponding prediction with wP^ [u> = 0) reasonably 
accounts for the width in the small-g regime. 

It is seen by comparing Figs. 1151 and 1181 that the AOP 
portion of the longitudinal current spectra for the smaller 
density p — 0.95 is located at lower frequencies and its in- 
tensity is more enhanced. These differences simply reflect 
the shift of the AOP to lower frequencies accompanied by 
an increase of its spectral intensity for p = 0.95. In addi- 
tion, the positive dispersion of the sound velocity and the 
deviation from the Lorenzian shape of the spectra start 



V. CONCLUDING REMARKS 

In this paper, we presented a theoretical investigation 
on the high-frequency collective dynamics in liquids and 
glasses at microscopic length and time scales based on 
the MCT for ideal liquid-glass transition. We focused 
on recently investigated issues from IXS and computer- 
simulation studies for dynamic structure factors and lon- 
gitudinal and transversal current spectra: the anomalous 
dispersion of the high-frequency sound velocity and the 
nature of the low-frequency excitation called the boson 
peak. The results of the theory are demonstrated for a 
model of the LJ system. 

The MCT explains the evolution of the structural- 
relaxation processes as precursor of the glass transition 
at a critical temperature T c , which is driven by the mu- 
tual blocking of a particle and its neighbors (cage ef- 
fect) |U|2^|. On the other hand, the theory also predicts 
the development of the harmonic vibrational excitations 
in stiff-glass states T < T c , called the AOP. The AOP 
is caused by the strong interaction between density fluc- 
tuations at microscopic length scales and the arrested 
glass structure, and exhibits the properties of the boson 
peak [24|. As demonstrated in Sec. IIII CI the AOP per- 
sists also near and above T c , and considerable portion of 
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the microscopic process in systems ranging from high-T 
liquids down to deep-in-glass states can be described as 
a superposition of harmonic vibrational dynamics. We 
investigated in Sec. lIIII how the interference of the sound 
mode with these structural-relaxation processes and vi- 
brational excitations shows up as anomalies in the sound- 
velocity dispersion. 

In the vicinity of T c , the structural-relaxation contri- 
butions to the sound-velocity dispersion show up at very 
small wave numbers, which can fully be appreciated only 
with the logarithmic q axis (cf. Fig. lll|) . Therefore, with 
the linear-g axis adopted in conventional studies on the 
dispersion relation, the sound velocity v q in the small- 
q limit is located above the hydrodynamic value Vq (cf. 
Fig-He)), and this even at T = 4.0 (cf. Fig.^b)) which 
is more than twice of the critical temperature T c . We no- 
tice that this temperature is also higher than the freezing 
temperature Tf » 3.3 which is estimated based on the so- 
called Hansen- Verlet criterion [22) • Thus, except for very 
high temperature like T — 10 (cf. Fig. Ufa)), the varia- 
tion of the sound velocity v q which is visible in the linear 
q axis basically reflects the contribution from the micro- 
scopic process. This theoretical result is in contrast to 
the traditional picture based on the viscoelastic model 
according to which the anomalous dispersion is fully as- 
cribed to the structural relaxation, but is in agreement 
with finding from recent IXS studies on simple liquid 
metals near the melting temperature [l^, LL-J OH OH • 

The description of the microscopic process in terms 
of the AOP - a superposition of harmonic oscillations 
of particles inside their cages - implies that the decay 
of the memory kernel at microscopic times log 10 t ^ — 1 
(cf. Fig. reflects the dephasing of different oscillatory 
components in the force fluctuations. This is in consis- 
tent with the implication from the computer-simulation 
study on harmonic glass Q. Such microscopic process, 
despite arrested structural relaxations, accounts for the 
positive dispersion in glass states (cf. Fig.[3J), whose pres- 
ence has been reported from recent MD and IXS stud- 
ies 0, S IS El fllf - That the microscopic process 
due to the AOP persists both in liquid and glass states 
with only weak temperature dependence of its character- 
istic time or frequency (cf. Fig. [HJ is in accord with the 
computer-simulation result presented in Ref. [9j. 

To see further implications of the theory, we show in 
Fig. [2] the sound velocities v q (T) as a function of T for 
fixed q = 10~ 2 (solid line) and q = 0.5 (dash-dotted 
line). (The argument T shall be added here to empha- 
size the T dependence.) The density is fixed to p = 1.093, 
and v q (T) are obtained from Eq. (|24b|l with A(u) de- 
termined for each temperature (cf. Fig. 0). Assuming 
that LJ particles under study are of argon size (« 0.34 
nm), the wave number q = 0.5 (1.5 nm" 1 ) is in the low- 
est momentum-transfer range in IXS measurements, and 
q = 10~ 2 (0.03 nm -1 ) is a typical momentum transfer in 
Brillouin light scattering (BLS) experiments. As shown 
in Fig.[3Jc) for a stiff-glass state T = 0.5, a negative dis- 
persion is observable for q < 1 where the sound velocity 
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FIG. 21: Sound velocities v q (T) as a function of T at 
p = 1.093 for q = 10~ 2 (solid line) and q = 0.5 (dash-dotted 
line), which are based on Eq. H24bt with A(u) determined for 
each temperature (cf. Fig.|SJ. Hydrodynamic sound velocities 
v (T) for T > T c and v (T) = v^(T) (see text) for T < T c are 
denoted as open circles and filled squares, respectively. v^(T) 
for T > T c , given by w£,(T) = «o( r )/\A ~ To, are also de- 
noted as filled squares. Dotted line shows the square-root sin- 
gularity translated to v£>(T), i.e., <,(T) = v (T)y/l - f (T) 
with f (T) = / C + hoy/a/(l-X) for T < T c and f (T) = f§ 
for T > T c . Dashed line exhibits the /3-relaxation contribution 
to Vq{T) which is based on Eq. I24bl with A(u>) determined 
from the MCT asymptotic formula I3LH . 



v q (T) is smaller than the hydrodynamic (q — * 0) value. 
This is due to the dominant harmonic nature of the dy- 
namics at T = 0.5, because of which there appears a 
frequency interval where A(a>) is smaller than its u> — > 
limit (cf. Fig. HOfl . The presence of the negative disper- 
sion is reflected in Fig. as the emergence of the low-T 
region where the IXS sound velocity (v q (T) for q = 0.5) 
is smaller than the BLS sound velocity (q = 10~ 2 ). Such 
T dependence of the IXS and BLS sound-velocity data 
in the low-T region has been reported for glycerol 0, 13^ 
(see Fig. 5(b) of Ref. fl). Thus, there is an experimen- 
tal result which is consistent with the presence of the 
negative dispersion predicted by the theory. That the 
difference between the IXS and BLS sound-velocity data 
in the low-T region is smaller in Fig. |2] than the one 
found in Fig. 5(b) of Ref. for glycerol and that such 
difference is not clearly observable in orthoterphenyl |39l | 
might be rationalized by the stronger influence of vibra- 
tional modes (boson peaks) in the dynamics of network 
glass former glycerol than in fragile glass formers like the 
LJ system and orthoterphenyl |40j. According to this 
reasoning, the presence of the negative dispersion is also 
expected in strong glass formers like silica where the in- 
fluence of vibrational dynamics is even stronger 40] . The 
sound velocity data for vitreous silica reported in Fig. 3 
of Ref. ^3 seems to indicate this possibility, but large er- 
ror bars do not allow ones to draw a decisive conclusion. 
It would therefore be of value to re-analyze the data for 
silica with improved resolution. 
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Hydrodynamic sound velocities vq(T) for T > T c and 
v (T) = v (T)/y/l - f (T) for T < T c are included in 
Fig. as open circles and filled squares, respectively. 
vq(T) increases with increasing T because in our study 
T is varied with density fixed, leading to increased pres- 
sures at elevated temperatures. On the other hand, the 
T dependence of vq(T) is dominated by that of fo(T) 
(cf. Fig.nj, and this explains the increase of io(T) with 
decreasing T. Here a side remark shall be added con- 
cerning the notion of the "hydrodynamic" sound veloc- 
ity v (T) for glass states referring to T < T c . The MCT 
in its idealized form, adopted in the present study, pre- 
dicts the structural arrest at the critical temperature T c 
which is located above the calorimetric glass-transition 
temperature T g . In reality there are slow dynamical pro- 
cesses - referred to as hopping processes - which restore 
ergodicity for T < T c . These processes, which are stud- 
ied within the extended version of the MCT |4l| > change 
the ideal elastic peaks nfq(T)5(ui) of density fluctuation 
spectra 4>g(u)) into quasielastic a peaks of nonzero width 
also for T < T c . f q (T) in real systems therefore has to 
be interpreted as an effective Debye- Waller factor |42| . 
which can be measured as area under the quasielastic 
a peak in or by determining the plateau of the 

4> q (t)-versus-log t curve. Correspondingly, Vo(T) should 
be interpreted as the extension of v^(T) - the sound 
velocity in the high-frequency limit of the a regime - in- 
troduced in Sec. llTTDl to T < T c . y ° (T) for T > T c are 
given by t£(T) - v (T)/^T=J$ (cf. Sec. HUE), and 
are also included with filled squares in Fig. [2] (T) 
defined below and above T c merge at T c . 

The mentioned dependence of u£,(T) on fo(T) has 
been utilized to extract experimentally the zero wave- 
number limit of the Debye- Waller factor via the rela- 
tion f (T) = 1 - [«o(T)/«° (T)] 2 HI E| and to test 
the square-root si ngularity as predicted by the MCT, 
fo(T) = / C + hoy/(T/(l-X) for T<T C and f Q (T) = / C 
for T >T C (cf. Eq. ©). The dotted line in Fig.l2"T1shows 
the square-root singularity translated to v^(T). In our 
system, there is no possibility to detect the singularity 
based on the IXS sound velocity (v q (T) for q = 0.5), 
which basically reflects the contribution from the mi- 
croscopic process (cf. Fig. 111(1 . On the other hand, 
the BLS sound velocity (q — 10~ 2 ) probes much lower 
momentum-transfer range, and is thus affected also by 
the structural-relaxation processes. In particular, the (3- 
relaxation contribution to v q (T) plays an important role 
in the vicinity of T c (cf. Fig. Illfl . as can also be inferred 
by comparing the solid and dashed lines in Fig. |^ the 
latter being obtained based on Eq. (|24b|) with A(ui) de- 
termined from the MCT asymptotic formula (|3L)[) for the 
(3 relaxation. Thus, the (3 relaxation must be included 
for reliable determination of fo(T) via the sound veloc- 
ity data, which is in full agreement with the conclusion 
drawn in Ref. 44]. (cf. Fig. 13 of Ref. where BLS 
data for CaKN03 corresponding to solid line and filled 
squares in Fig. 1211 are shown.) 

We studied in Sec. IIVI how the AOP manifests itself 



in the spectral shape of the glass-state longitudinal and 
transversal current spectra at low frequencies. For small 
wave numbers where the sound resonance frequencies are 
located below the low-frequency threshold of the AOP, 
^L(T)max ^ ^ ^ there is practically no effect from the 
AOP, and the longitudinal and transversal current spec- 
tra are dominated by the sound excitations of the Loren- 
zian spectral shape. By increasing the wave numbers 
so that the sound resonance frequencies enter the AOP 
region, f2_ < w^ T ' )max , the hybridization of the sound 
mode with the AOP becomes important, and the shape 
of the current spectra deviates from Lorenzian. As the 
wave number is increased further, the interference of the 
sound mode with the AOP leads to the positive disper- 
sion of the sound velocities, pushing ujq max consider- 
ably above frequencies expected from the bare disper- 
sion. Thereby, there opens a frequency window, where 
the appearance of the AOP itself can be observed in the 
spectra at low frequencies. Furthermore, our approxi- 
mate formulas (I36H predict that (i) there is no q depen- 
dence in the peak frequency waop of the AOP portions 
in the current spectra, (ii) waop is nearly the same for 
both of the longitudinal and transversal current spectra 
since rh"(u>) ~ m^(u>) as shown in Fig.0^c), and (iii) the 
intensity of the AOP portion is stronger for the transver- 
sal current spectra than for the longitudinal ones in the 
pseudo first Brillouin zone. These theoretical predictions, 
demonstrated in Figs. to |201 for the LJ system, are 
in agreement with previous findings from MD and IXS 
studies [1, 0, 0, 0| . Thus, the low-frequency excita- 
tions observable at the same resonance frequency u>aop 
in both of the glass-state longitudinal and transversal 
current spectra are the manifestation of the AOP. Usu- 
ally, the boson peaks are discussed below T g . On the 
other hand, the glass state in the present paper refers 
to T < T c . Since T c is located above T g , our theoretical 
results for glasses apply also for some ranges of the liquid 
state in the conventional terminology. 

According to our theoretical results, it is essential to 
have the positive dispersion of the transversal sound ve- 
locity for the appearance of the low- frequency excitations 
in the transversal current spectra, whose presence has 
never been discussed so far. This comes out from our 
study on the density dependence of the dispersion re- 
lation, shown in Figs. I17f a) and I20f a^l for the higher 
and the lower densities, respectively. Such density de- 
pendence of the dispersion relation is in accord with the 
one reported from the computer simulation for vitreous 
silica (cf. Fig. 4 of Ref. It would therefore be of 

value if simulation studies could test whether or not our 
theoretical prediction agrees with their data. 

We notice that the intensity of the low-frequency exci- 
tations for the LJ system studied in this paper is weaker 
than the one, e.g^ known from computer-simulation re- 
sults on silica |sl.ll0| : the low- frequency excitations in the 
present study are shoulders rather than peaks. Again, 
this feature may be ascribed to the fact that the LJ sys- 
tem can be classified as a fragile glass former, in which 
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FIG. 22: (a) The memory-kernel spectrum rh"(fj) and (b) 
the SGA spectrum i?i5 GA (u) obtained with the cutoff q* — 40 
(solid lines), 80 (dashed lines), and 120 (dotted lines) for 
the hard-sphere system at packing fraction ip = 0.6. As 



ated within the Percus-Yevick approximation, and the units 
of length and time have been chosen so that the hard-sphere 
diameter d = 1 and the thermal velocity v = 2.5. 



the boson peak is less pronounced than in strong glass 
formers like silica [4p|. 

Finally, let us make a comment on the cutoff prob- 
lem mentioned in Ref. [24j. As already noticed there, 
the MCT results for the hard-sphere system in stiff-glass 
states change if the cutoff q* of the wave-number grids 
used in the calculations is varied. This is demonstrated 
in Fig. I22f a) showing the q* dependence of the memory- 
kernel spectrum m"(w) for the hard-sphere system at 
packing fraction ip = 0.6 studied in Ref. |2J|. It is 
seen that the intensity of m"(u>) increases with increasing 
q* . Furthermore, the broad harmonic-oscillation "peak" 



reflecting the AOP, observable in m"(w) for the cutoff 
q* = 40 adopted in Ref. [24j], disappears for q* = 120; it 
shows up only as a shoulder. Such q* dependence reflects 
the slow decrease towards zero of the direct correlation 
function Cfc of the hard-sphere system for k tending to 
infinity, because of which the relevant coupling coeffi- 
cient 14 (cf. Eq. (|2*c|) ) does not approach zero even in 
the k — > oo limit. However, this does not mean that 
there is no AOP in the hard-sphere system. As discussed 
in Sec. 1111 CI the AOP portion of m"(ui), even when it 
is buried under the tail of the quasielastic peak, is well 
extracted by the SGA spectrum mg GA (w). Figure 122T b) 
demonstrates that there is no q* dependence in Wg GA (w) 
provided a sufficiently large q* is chosen. In this sense, 
the AOP is well defined in the hard-sphere system. But, 
it does not show up as a peak in m"(uj) due to the two- 
mode contributions, which are neglected within the SGA. 
Thus, the increase of the intensity of rh"(ui), in particu- 
lar, the smearing out of the "gap" for < u) < f2_, with 
increasing q* is due to the development of the two-mode 
contributions. Because of the strong anharmonic effects, 
no harmonic-oscillation peak in the spectrum m"(u>) is 
expected for the hard-sphere system. 

On the other hand, introducing the cutoff is equivalent 
to softening the hard-sphere potential. Indeed, we have 
confirmed that there is no cutoff problem in the LJ sys- 
tem considered in this paper, and that all the essential 
results presented in Ref. [24( can be reproduced for this 
system. In particular, one observes harmonic-oscillation 
peak for the LJ system in stiff-glass states as shown in 
Fig. [7Ja) , which is free from the mentioned cutoff prob- 
lem. Thus, the whole physics discussed in Ref. [24( for 
stiff-glass states remain valid for systems where particles 
interact with regular interaction potentials. 
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